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1.0  INTRODUCTION 


1.1  Background 

Advanced  antiair  missile  guidance  systems  will  make 
extensive  use  of  infrared  (IR)  sensors  to  counter  ra¬ 
dar  cross-section  reduction  and  standoff  jamming,  to 
achieve  lower  miss  distance,  and  to  increase  high- 
altitude  missile  dynamic  performance.  IR  sensors  will 
also  provide  an  added  discrimination  capability  to  de¬ 
feat  self-screening  countermeasures  and  select  targets 
from  closely  spaced  formations.  The  Wide-Area  Guid¬ 
ance  and  Control  Program  is  tasked  with  developing 
guidance  technology  for  use  in  advanced,  long-range, 
surface-to-air  missiles  (SAM).  Such  missiles  must  be 
“self-targeting,”  that  is,  the  missile  guidance  conducts 
wide-field  angle  search  to  compensate  for  poor  tar¬ 
geting  (i.e.,  midcourse  flight)  accuracy  against  long- 
range  targets.  While  advances  in  detector  technology 
potentially  offer  very  sensitive  IR  seeker  systems,  de¬ 
tection  performance  will  most  likely  be  limited  by  the 
background  clutter. 

An  investigation  of  clutter-imposed  limitations  to 
missile  IR  seeker  performance  is  being  conducted  to 
determine  the  likely  performance  of  a  guidance  sen¬ 
sor  that  could  be  used  on  long-range  SAMs.  The  ob¬ 
jective  of  this  effort  is  to  develop  clutter  processing 
algorithms  and  to  apply  those  processing  techniques 
to  data  in  order  to  determine  the  clutter-limited  per¬ 
formance  of  advanced  IR  seekers.  The  principal  clut¬ 
ter  processing  technique  investigated  here  is  a 
combination  of  digital  filtering  followed  by  a  pulse 
discrimination  algorithm. 

1.2  Missile  Infrared  Seekers 

A  long-range  SAM  with  IR  guidance  has  several 
unique  requirements  that  dictate  the  choice  of  signal 
processing.  First,  an  IR  seeker  must  acquire  and  track 
the  target  autonomously.  Potential  false  alarms,  if  not 
immediately  rejected,  result  in  the  waste  of  an  expen¬ 
sive  missile.  Second,  an  IR  seeker  must  be  capable  of 
rapid,  wide-field  search  to  compensate  for  significant 
targeting  and  midcourse  guidance  errors  that  may  oc¬ 
cur  when  targets  are  engaged  at  long  range.  Third,  the 
long-range  nature  of  the  missile  requires  sufficient 
seeker  acquisition  range  to  acquire  targets  from  the 
optimum  flight  altitude  to  allow  time  for  descent  to 
target  altitude.  Hence,  the  capability  to  detect  small 
signature  targets  in  terrain  and  cloud  background  clut¬ 
ter  at  long  range  is  crucial  for  mission  success. 

Achievement  of  high  IR  seeker  performance  will  be¬ 
come  significantly  more  important  as  enemy  capabil¬ 


ities  evolve.  Autonomous  IR  seeker  acquisition 
becomes  more  important  as  target  radar  cross  section 
is  reduced  and  the  IR  mode  assumes  greater  initial  de¬ 
tection  responsibility.  Initially,  the  IR  seeker  will  have 
primary  acquisition  responsibility  against  high-speed, 
higher  altitude  (i.e.,  higher  IR  signature)  targets  with 
the  missile  seeker  radar  having  the  primary  responsi¬ 
bility  for  detecting  lower  speed  and  lower  altitude  (i.e., 
lower  IR  signature)  targets.  Ultimately,  an  IR  seeker 
capability  against  low-signature  targets  is  desired. 

Because  of  the  severity  of  the  targeting  and  mid¬ 
course  guidance  requirements  for  a  long-range  mis¬ 
sion,  the  enemy  will  develop  tactics  and  electronic 
countermeasures  (ECM)  to  degrade  these  functions. 
As  targeting  and  midcourse  guidance  quality  is 
reduced,  the  search  capability  of  all  missile  guidance 
modes  must  be  increased.  Increases  in  the  size  of  the 
search  field  of  regard,  the  search  angular  scan  rate, 
and  acquisition  range  are  needed  to  offset  the  effects 
of  poorer  target  localization.  The  long-range  missile 
must  have  acquisition  and  intercept  capability  even  in 
situations  when  target  range  is  unknown  or  external 
communications  (for  midcourse  guidance)  are  lost. 
Uncertainties  in  range  also  require  longer  search  time, 
hence  a  lower  false-alarm  rate. 

As  external  targeting  becomes  poorer,  particularly 
when  target  range  information  is  degraded  or  unavail¬ 
able,  the  need  for  the  missile  to  fly  a  propulsion- 
efficient  trajectory  increases.  Since  propulsion-efficient 
trajectories  result  in  high-altitude  flight,  range  to 
potential  targets  is  increased.  The  increased  range  re¬ 
quirement,  coupled  with  the  need  for  more  search 
capability  and  the  need  to  detect  lower  signature  tar¬ 
gets,  results  in  an  IR  seeker  requirement  for  signifi¬ 
cantly  greater  sensitivity. 

Advances  in  detector  technology  offer  large  in¬ 
creases  in  sensor  sensitivity,  e.g.,  two  orders  of  mag¬ 
nitude  over  current  systems.  Increased  IR  sensitivity 
offers  the  potential  for  the  needed  long-range  detec¬ 
tion  of  low-signature  targets.  However,  actual  capa¬ 
bility  depends  on  the  development  of  satisfactory 
techniques  for  clutter  suppression. 

In  addition  to  the  unique  long-range  mission  require¬ 
ments  discussed  above,  all  IR  seekers  have  common 
constraints  that  affect  their  clutter  discrimination  capa¬ 
bility.  Most  important  are  resolution,  stabilization,  and 
high  closing  velocity.  Typically,  IR  seekers  have  an 
aperture  of  about  5  cm  resulting  in  a  typical  minimum 
resolution  of  0.5  mrad  at  5  fim  and  1.0  mrad  at  10 
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nm.  Current  IR  seekers  have  a  typical  stabilization 
drift  rate  of  ~0.1  degree/sec. 

There  are  three  basic  techniques  for  point-target  dis¬ 
crimination  in  clutter.  The  first  is  the  use  of  a  reticle 
to  encode  point-target  data  in  a  background  scene. 
This  approach  allows  continual  observation  of  a  wide 
field  with  a  single  detector  while  a  spatially  small  tar¬ 
get  is  tracked.  The  technique  is  used  in  most  current 
IR  seekers.  Unfortunately,  the  reticle  modulates  any 
background  structure  within  the  field  of  view  (FOV). 
This  unwanted  modulation  can  interfere  with  target 
detection,  particularly  if  greater  seeker  sensitivity  is 
required.  An  even  greater  drawback  to  reticle-based 
systems  is  the  great  difficulty  in  separating  multiple 
targets  or  targets  and  decoys. 

The  drawbacks  of  the  first  technique  are  mitigated 
by  the  second  technique,  i.e.,  scanning  a  larger  FOV 
with  small  instantaneous  field-of-view  (IFOV)  dei<.c- 
tors.  Scanning  systems  can  be  implemented  with  sin¬ 
gle  detectors  or  multiple  detectors  up  to  large  arrays. 
Detector  size  is  matched  to  optical  resolution  (i.e.,  to 
the  size  of  the  blur  circle).  Discrimination  is  based  on 
the  fact  that  most  of  the  scene  power  is  concentrated 
in  large  features.  Small  features,  on  the  scale  of  the 
blur  circle  size,  remain  a  source  of  clutter.  As  resolu¬ 
tion  increases  (blur  circle  and  detector  size  decrease), 
the  clutter  power  entering  the  system  is  reduced.  Flence 
detection  capability  is  improved  as  resolution  increases. 

The  third  technique  is  the  staring  system.  Again, 
small  IFOV  detectors  are  used.  Discrimination  is  based 
on  moving  target  detection.  Each  detector  looks 
(stares)  at  a  portion  of  the  background.  Successive 
frames  are  subtracted  to  detect  the  motion  of  a  target 
from  one  pixel  to  the  next.  Ideally,  all  background 
power  is  removed  in  the  subtraction  process,  leaving 
a  result  that  is  corrupted  only  by  system  noise.  Un¬ 
fortunately,  there  are  practical  implementation  prob¬ 
lems,  including  changes  in  background  power  during 
the  stare  period,  caused  by  missile  motion  relative  to 
the  background  (e.g.,  projected  velocity  on  the  three¬ 


dimensional  scene),  by  seeker  stabilization  (drift)  er¬ 
rors  resulting  in  pixel  misregistration,  and  by  temporal 
variations  in  the  scene.  The  relatively  large  IFOV  of 
a  missile  seeker  (compared  to  advanced  search  set)  re¬ 
quires  a  longer  stare  time  for  detection.  The  longer 
stare  time,  coupled  with  background  changes,  drift, 
and  the  need  to  scan  a  wide  search  field  rapidly,  makes 
implementation  of  staring  processing  unlikely  in  a 
long-range  missile  IR  seeker. 

For  the  reasons  shown,  this  report  concentrates  on 
development  and  test  of  clutter  suppression  algorithms 
for  scanning  IR  seekers. 

1.3  Approach  to  Algorithm  Development 

The  background  data  used  in  this  study  are 
described  first.  A  characterization  of  the  background 
imagery  is  then  provided.  The  processing  techniques 
used  in  the  analysis  are  then  described,  followed  by 
the  results  of  processing  the  data  by  these  techniques. 

The  description  of  background  data  includes  a  sum¬ 
mary  of  data  requirements  and  a  survey  of  available 
IR  clutter  data.  Reasons  for  selecting  the  background 
imagery  used  in  this  study  and  a  detailed  description 
of  that  data  are  provided. 

Background  data  are  then  characterized  and  ar¬ 
tifacts  in  the  data  examined  and  corrected.  The  char¬ 
acterization  is  used  to  guide  the  implementation  of 
signal  processing  algorithms.  In  some  cases,  the  back¬ 
ground  characterization  can  be  used  to  estimate  seeker 
performance  directly  and  to  make  a  comparison  be¬ 
tween  scenes  and  spectral  bands. 

Finally,  the  original  background  data  are  directly 
processed  with  clutter-suppression  algorithms  to  de¬ 
termine  the  effectiveness  of  the  algorithm  and  ultimate 
clutter-limited  performance  of  an  advanced  seeker. 
These  algorithms  are  described  and  the  reasons  for 
their  selection  explained.  Comparisons  of  sensor  per¬ 
formance  are  made  between  different  spectral  bands 
and  different  scene  types.  Conclusions  are  in  the  form 
of  expected  sensor  performam.^ 
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2.0  BACKGROUND  DATA 


2.1  Data  Requirements 

Background  data  needed  to  support  a  wide-area 
missile  IR  seeker  are  desired  in  the  form  of  high  spa¬ 
tial  resolution,  two-dimensional,  contiguous  back¬ 
ground  maps  (scenes).  Since  the  data  support  the 
evaluation  of  a  scanning  system,  high  relative  radio- 
metric  accuracy  along  one  of  the  scene  dimensions 
(scan  direction)  is  needed.  Oversampled  data  in  this 
direction  are  needed  so  filtering  can  be  applied.  Ab¬ 
solute  radiometric  accuracy,  to  within  20%,  is  required 
to  provide  a  good  estimate  of  sensor  performance. 
Scanning  background  measurement  sensors  are  most 
likely  nece,ssary  to  provide  the  radiometric  accuracy 
and  oversampling  needed  to  predict  performance  of 
a  scanning  system. 

The  wide-area  missile  mission  requires  scenes  with 
down-looking  geometries.  Typical  seeker  search  is 
from  high  altitude,  looking  below  the  horizontal  at 
depression  angles  of  5  to  20°.  Since  a  representative 
set  of  all  scenes  of  interest  is  an  enormous  quantity 
of  data,  this  study  concentrates  on  cloud  scenes. 
Clouds  present  a  variety  of  difficult  clutter  conditions, 
particularly  over  water. 

Background  data  in  two  spectral  bands  are  need¬ 
ed:  middle-wavelength  data  and  long-wavelength  data. 
The  middle-wavelength  data,  nominally  4  to  5  ^jm,  re¬ 
quire  an  angular  spatial  resolution  of  approximately 
0.5  mrad.  Long-wavelength  data,  nominally  8  to  12 
^im,  require  an  angular  spatial  resolution  of  about  1 .0 
mrad.  Data  taken  at  altitudes  lower  than  typical  wide- 
area  missile  search  altitude  (60,000  to  100,000  ft)  can 
be  obtained  with  lower  resolution  instruments,  i.e., 
instruments  that  have  the  equivalent  linear  spatial  reso¬ 
lution.  Acquisition  of  simultaneous  middle-  and  long- 
wavelength  data  is  preferred  to  make  a  sidc-by-sidc 
comparison  of  performance  between  these  two  ma¬ 
jor  spectral  hands. 

The  noise  in  the  background  measurement  sensor 
must  be  below  that  expected  for  advanced  sensors. 
Sensor  radiometric  accuracy  of  10  W/em'/sr 
(middle  wavelength)  and  10  W/em’/sr  (long  wave¬ 
length)  is  needed.  The  data  resolution  (digital  sample 
increment)  must  be  less  than  or  equal  to  the  rms  noise 
requirements  so  that  the  sensor  noise  rather  than 
digitizing  noise  is  the  predominant  noise  source. 

Each  scene  should  have  sufficient  samples  to  pro¬ 
vide  an  estimate  of  false-alarm  rate.  A  minimum  of 
lO'’  pixels  is  necessary;  lO**  pixels  or  more  arc  desired. 


Table  1  summarizes  the  requirements  for  back¬ 
ground  data.  These  requirements,  while  designed  to 
support  development  of  a  scanning  IR  seeker  for  a 
wide-area  missile,  are  believed  to  be  equally  applica¬ 
ble  to  scanning  infrared  surveillance  from  aircraft. 


2.2  Data  Survey 

A  survey  of  available  IR  background  data  was  con¬ 
ducted  to  gain  an  understanding  of  the  available  back¬ 
ground  data  and  to  guide  selection  of  the  data  for  use 
in  this  study.  Table  2  summarizes  some  of  the  avail¬ 
able  data  sources.  Comments  on  these  sources  are  giv¬ 
en  below. 

The  Armament  Laboratory  at  Eglin  Air  Force  Base 
maintains  two  databases:  BASES  and  TABILS.  The 
BASES  system  is  designed  to  support  the  development 
of  missile  seekers.  It  consists  of  a  pod  mounted  on 
an  F4  aircraft  with  3  to  5  and  8  to  I2fim  imagery  (1.0 
mrad)  capabilities.  These  data  are  collected  for  a 
specific  program,  or  SPO,  and  are  not  well  cataloged. 
Targets  are  primarily  tactical  aircraft  (all  aspects)  at 
close  range.  The  background  data,  at  ranges  of  in¬ 
terest,  do  not  have  sufficient  signal-to-noise  ratio  for 
use  in  this  study. 

The  .second  source  is  the  Target  and  Background 
Information  Library  System  (TABILS),  which  is  ac¬ 
cessible  through  remote  computer  terminals  to  DoD 
agencies  and  contractors.  This  database  was  developed 
to  support  air-to-ground  tactical  battlefield  programs 
such  as  Assault  Breaker,  Tank  Breaker,  etc.  A  user’s 
manual  is  available.' 

The  Environmental  Research  Institute  of  Michigan 
(ERIM)  has  collected  a  large  volume  of  airborne  im¬ 
agery  with  their  M5  and  M7  line  scanners  as  part  of 
the  Navy  Optical  Signatures  Program  as  well  as  for 
other  sponsors.  This  work  started  in  1976  with  an  ex¬ 
tensive  survey  of  available  background  data,"  fol¬ 
lowed  by  data  collection  of  ground  scenes  from  low 
altitude.  Typical  data  arc  collected  at  1000  to  2000  ft 
altitude  in  three  IR  spectral  bands.  Detector  resolu¬ 
tion  of  2  to  2.5  mrad  gives  high  spatial  resolution  on 
the  ground.  .A  summary  of  the  data  collected  for 


Air  I  orLC  Armameiil  I  aboraiorv ,  rurvi’i  and  Hack^roiind  In- 
iDniHition  I  dtrury  Sv-fcni  {  rAlil/  S]  L  sit's  Munuul  {Way  I. 
IWS:). 

.1.  Hoard.  ,1  Brailuailc,  and  R.  7  iiriior.  Infrared  Buck  ground 
Survey  and  Analysis.  I  RIM  1IS0()0-1-I  (jim  1976). 
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Table  1  Requirements  for  wide-area  missiie  IR  background  data. 


Parameter 

Requirement 

Geometry 

60,000-l(X),000  ft  altitude 

5-20*  depression  angle 

Scenes 

Cloud: 
solid  layers 
multiple  layers 
broken  clouds 
storms  and  thunderheads 

Land 

Spectral  bands 

4.0-5 .0  )jm 

8.0-12  fim 

Angular  resolution 

0.5  mrad 

1 .0  mrad 

Sensitivity 

10“^  W/cmVsr 

10'*  W/cmVsr 

Radiometric  dynamic  range 
(from  zero) 

12  bits 

12  bits 

Radiometric  accuracy 

better  than  20®7o 

Oversampling 
(scan  direction) 

2:1  or  greater 

Number  of  pixels  per  scene 

>  10* 

1978-1979,  as  well  as  earlier  data,  is  given  in  Ref.  3. 
More  recently  (1983),  long-wavelength  European  ter¬ 
rain  background  data  were  collected. 

A  number  of  government  agencies  and  contractors 
are  tasked  to  collect  data  for  the  long-range  detection 
of  air  targets.  General  Electric  and  ITT,  under  con¬ 
tract  to  the  Air  Force  Weapons  Laboratory,  have  de¬ 
veloped  middle-wavelength  airborne  IR  search  and 
track  systems.  These  systems  were  tested  at  Eglin  Air 
Force  Base  in  the  summer  of  1985.  Data  produced  by 
the  tests  were  not  available  for  this  analysis.  The  Navy 
is  planning  to  continue  the  development  with  full-scale 
engineering,  starting  in  late  1986. 

The  Navy  Background  Measurements  and  Analy¬ 
sis  Program  (BMAP)  is  being  conducted  as  a  succes- 


'  J.  LaRocca  and  D.  J.  Witte,  Handbook  of  the  Statistics  of 
Various  Terrain  and  Water  (la)  Backgrounds  from  Selected  U.S. 
Locations.  ERIM  139900-1-X  (Jan  1980). 


sor  to  the  background  portion  of  the  Optical 
Signatures  Program.  The  program  is  tasked  with  col¬ 
lecting  data  needed  to  design  surface-  and  air-based 
IR  search  and  track  systems.  The  initial  instrument 
is  a  high-spatial-resolution,  high-sensitivity  device  with 
simultaneous  4  to  5  /tm  and  8  to  1 1  ^m  coverage. 
Initial  data  collection  began  with  surface  measure¬ 
ments  in  1984  (Ref.  6)  and  has  continued  with  airborne 
(P3  aircraft)  measurements.  A  follow-on  instrument 


■•r.  L.  Lucke,  el  al.,  “The  Navy’s  Infrared  Background  Mea¬ 
surement  and  Analysis  Program,”  in  Proc.  IRIS  Specialty  Group 
Meeting  on  Targets,  Background,  and  Discrimination,  San  Diego 
(Feb  11-14,  1986). 

'R.  a.  Steinberg,  “Navy  IR  Background  Measurements  and 
Analysis  Program,”  in  Proc.  Tri-Service  Infrared  Backgrounds 
Symp.,  Oct  18-20,  1983,  Mitre  Corporation,  pp.  216-223. 

'^A.  Hirschman,  “BMAP  Surface-Based  Background  Measure¬ 
ment  Activities  at  Montauk  Point,  New  York,  August  1983,” 
Naval  Surface  Weapons  Center,  White  Oak,  Silver  Spring,  Md., 
draft  report  (Nov  1984). 
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Table  2  Summary  of  available  background  data. 


Source 

agency 

Database 

Platform 

Sensors 

AFATL 

BASES 

TABILS 

F4 

Various 

AGA  3-5,  8-12  /tm 
Imagery 

AFFTD 

KC-135 

Spectrometer,  3-5  /jm 
FLIR,  8-12  ^m 

AFWAL 

FLIR,  8-12  tim 

AFGL 

NKC-135 

Spectrometer,  3-5  #im 
NORELCO,  3-5 

FLIR,  8-12  Mm 

NSWC/DL 

Ground- 

based 

Spectrometer,  3-5  pm 
FLIR,  3-5  and  8-12  Mm 

NRL 

BMAP 

P-3 

Radiometer,  3-5  and 

8-12  pm 

NWC 

Target 

Signatures 

Program 

Various 

DARPA/Lockheed 
(ERIM  library) 

BMP,  CAMP 
HICAMP 

U2 

Scanner;  several  staring 
imagers,  10-12  Mm 

General 

Dynamics 

Ground- 

based 

Spectrometer,  3-5  and 
10-12  Mm 

Block 

Engineering 

Ground- 

based 

Spectrometer,  3-5  pm 

General  Electric 
Syracuse 

C-54 

Imager,  8-12  Mm 

ERIM 

C-47 

M-5  and  M-7  multi- 
spectral  scanner 

MIT/LL 

FLIR,  3-5  and  8-17  Mm 

JHU/APL 

P-3 

Radiometer,  3-4  and 

10-12  Mm 

with  an  expanded  number  of  detectors  per  spectral 
band  (from  16  to  120)  and  higher  spatial  resolution 
(0.25  mrad)  is  under  development. 

General  Dynamics,  NSWC/DL,  and  Block  En¬ 
gineering  have  ground-based  measurement  systems. 
These  organizations  maintain  excellent  calibrated  sig¬ 
nature  libraries  of  tactical  aircraft  flown  at  low  alti¬ 
tudes.  A  summary  of  spectral  background  measure¬ 


ments  made  by  General  Dynamics  is  given  in  Ref.  7, 
and  a  compendium  of  their  target  signatures  library 
is  given  in  Ref.  8. 


(ieneral  DynamitS' Pomona  Division,  Background  Spectral 
Radiance  and  Contrast  in  the  \ear-i'\',  Mid-IR.  and  I  H  IR 
Regions.  TM6-I25PH-4.t8  (Ocl  IV76). 

"Cieneral  Dynamics/Pomona  Division,  Availahle  Target  Infrared 
Spectral  Signature  Data.  TM6-I25PH-4SO  (Sep  1981). 
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The  Applied  Physics  Laboratory  has  measured  the 
radiometric  properties  of  sea  surface  with  low-altitude 
(SOO  to  3000  ft)  flights  of  a  line-scanning,  dual-band 
(3.6  to  4.1  nm  and  10.1  to  11.1  ;tm  radiometer.  This 
instrument  has  high  sensitivity  (NEAT  of  0.003  K  in 
each  band)  and  low  spatial  resolution  (10  to  2(X)  mrad). 
An  instrument  description  and  typical  data  are  given 
in  Ref.  9. 

Another  major  source  of  target  and  background 
data  is  the  DARPA  Infrared  Data  Library  located  at 
ERIM.  This  library  serves  as  the  archival  collection 
of  data  with  secret  security  classification.  Major  pro¬ 
grams  represented  in  the  library  are  the  DARPA- 
sponsored  Background  Measurements  Program 
(BMP),  the  Calibrated  Airborne  Measurements  Pro¬ 
gram  (CAMP),  and  the  Highly-Instrumented, 
Calibrated  Airborne  Measurements  Program  (both 
HICAMP  and  HICAMP-II).  A  report  on  library 
holdings'"  lists  the  contents  of  the  library. 

The  Air  Force  Geophysics  Laboratory  (AFGL) 
operates  an  instrumented  NKC-135  aircraft  to  acquire 
target  and  background  signature  data.  Over  the  past 
10  years  AFGL  has  collected  in-flight  spectral  and  spa¬ 
tial  data  on  several  tactical  and  strategic  aircraft  (e.g., 
T-38,  F4,  F5,  F14,  F15,  A4,  A7,  Bl,  B52,  C140, 
KC135,  etc.).  The  measurement  aircraft  typically  flies 
at  6  km  (20,0(X)  ft)  altitude,  but  is  capable  of  acquir¬ 
ing  data  at  altitudes  ranging  from  0.3  to  13  km  (1000 
to  42,000  ft).  Data  can  be  taken  in  the  down-,  side-, 
forward-,  and  rear-looking  modes. 

This  aircraft  contains  a  number  of  IR  measurement 
instruments  for  various  purposes.  Spectral  measure¬ 
ments  are  made  with  Michelson  interferometers  hav¬ 
ing  spectral  resolution  of  0.1  (for  plume)  to  32  cm  “ 
Imagery  data  are  taken  in  two  broad  spectral  bands 
(3  to  5,  8  to  12  ^m)  and  subbands  within  these  inter¬ 
vals.  The  mid-wave  IR  (MWIR)  imagery  data  are  taken 
at  1  mrad  resolution  using  an  AC -coupled  system. 
Consequently,  the  data  only  provide  differential  radi¬ 
ance  values.  The  8  to  12  fim  forward-looking  IR 
(FLIR)  system  consists  of  two  9-element  HgCdTe  de¬ 
tectors  with  serial  scan  at  TV-compatible  rates. "  The 
IFOV  is  approximately  0.5  mrad,  and  a  typical  im¬ 
age  consists  of  192  x  256  pixels.  Individual  frames 
can  be  spliced  together  to  form  a  long-track  mosaic 
on  the  order  of  0.75  km  by  40  km.  Two  such  tracks 
have  been  constructed:  one  in  the  American  West  (near 


‘'K.  Peacock.  Radiometric  Temperature  Characteristics  of  the 
(Xean  Surface,  JHIJ/APL  STD-R-.t77  (Oct  1980). 

'"A.  J.  l.aRocca  and  B.  A.  Wtialev,  DARPA  Infrared  Data  Li¬ 
brary.  ERIM  l55300-24-.\  (Jul  i984). 

"J.  Shumsky  and  J.  H.  Schummcrs,  AFGL  El. IR  NKC-/.t5  Air¬ 
craft  Instrument.  AFGL  TR-82-0031  (Jan  26.  1982). 


Santa  Cruz,  California)  and  the  other  in  northern- 
latitude  ice  fields. 

The  AFGL  has  not  been  specifically  tasked  to  ac¬ 
quire  data  for  the  IR  search  and  track  (IRST)  (long 
range  detection)  problem.  However  they  have  acquired 
data  relevant  to  this  effort  in  the  course  of  scheduled 
missions.  These  include  horizontal  viewing,  2°  above 
and  below  the  horizon,  cloud  orbits,  and  cloud  fly- 
throughs.  These  data  have  not  been  cataloged. 

2.3  Selected  Data 

Selected  data  were  based  on  criteria  of  (a)  availabil¬ 
ity,  (b)  measurement  with  a  scanning  instrument,  and 
(c)  no  classification  restrictions.  Once  algorithms  are 
developed,  they  will  be  tested  on  classified  as  well  as 
unclassified  data. 

Three  data  sets  are  used  in  this  study:  images  from 
AFGL,  ERIM,  and  BMAP.  AFGL  scenes  are  all  long- 
wavelength  scenes  of  clouds  or  clouds  and  terrain 
viewed  from  high  altitude.  ERIM  scenes  are  short-, 
middle-,  and  long-wavelength  scenes  of  many  types 
of  terrain  as  well  as  of  water  taken  from  low  altitude. 
BMAP  data  consist  of  middle-  and  long-wavelength 
scenes  taken  from  the  ground.  Airborne  BMAP  data 
have  been  taken  and  are  expected  to  be  used  for  fu¬ 
ture  work.  Each  of  these  data  sets  is  described  in  the 
following  paragraphs. 

Two  sets  of  BMAP  data  were  supplied  by  the  Na¬ 
val  Research  Laboratory.  The  BMAP  instrument  con¬ 
sists  of  16  middle- wavelength  and  16  long- wavelength 
detectors  each  with  a  0.33  by  0.33  mrad  IFOV.  Data 
in  the  scan  direction  are  sampled  every  0.096  mrad 
resulting  in  an  oversampling  of  3.3  to  1 .  The  data  are 
stored  as  16-bit  words.  To  convert  the  data  to  radi¬ 
ance,  the  data  were  multiplied  by  the  multiplicative 
factor  (MF)  corresponding  to  the  correct  IR  band.  (For 
the  BMAP  data,  the  MF  equaled  the  noise  equivalent 
irradiance  (NEI)  of  the  system.)  Appendix  A  contains 
more  detailed  information  on  the  BMAP  instrument 
and  imagery  including  the  spectral  response  and  the 
impulse  response  of  the  sensor. 

The  BMAP  sensor  has  two  modes  of  operation. 
Data  can  be  collected  in  a  format  of  8  lines  of  long- 
wavelength  IR  followed  by  8  middle-wavelength  IR 
lines  per  frame.  In  this  8/8  format,  the  data  are  ex¬ 
tracted  from  coincident  detectors.  The  second  format 
is  16  lines  of  either  middle-  or  long-wavelength  IR  data 
per  frame.  In  this  mode,  the  frames  usually  alternate 
between  the  spectral  bands. 

APL  has  received  data  in  both  of  the  formats.  Fif¬ 
ty  frames  of  data  were  collected  at  Montauk  Point, 
Long  Island  in  the  8  long-wavelength  8  middle- 
wavelength  format.  The  background  scene  for  these 
data  was  overcast  sky. 
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The  second  set  of  BMAP  data  received  by  APL  was 
collected  at  Bedford,  Massachusetts.  These  data  were 
in  the  format  of  16  lines  of  alternating  spectral  bands 
per  frame.  APL  used  the  first  61  frames  of  these  data. 
Thirty-one  frames  were  mid  IR  and  30  frames  were 
long  IR.  Long-wavelength  IR  frames  only  had  IS  lines 
of  data  per  frame.  A  hole  in  a  cloud  deck  was  the  back¬ 
ground  scene  for  these  frames.  More  detail  on  the  im¬ 
ages  is  included  in  Appendix  A. 

Data  from  AFGL  were  the  second  set  used  in  this 
study.  The  instrument  is  a  two-detector  array,  serial 
scan,  FLIR  unit  with  a  10-position  filter  wheel  that 
can  select  various  spectral  bands.  All  data  supplied 
to  APL  were  taken  with  filter  I9005A  (7.84  to  13.16 
/tm).  The  FOV  on  this  seeker  was  0.36  (H)  x  0.43 
(V)  mrad,  but  the  digitizing  sample  size  was  0.S2  (H) 
by  0.45  (V)  mrad.  Eight  of  the  frames  have  an  addi¬ 
tional  telescope  decreasing  the  IFOV  to  less  than  0.1 
mrad.  The  sample  size  means  that  these  images  were 
slightly  undersampled  (as  opposed  to  the  BMAP  im¬ 
ages).  Video  data  (El A  RS-170)  are  recorded  and  sub¬ 
sequently  digitized  before  distribution  to  APL.  A  more 
detailed  description  of  the  AFGL  data  is  included  in 
Appendix  B. 

The  AFGL  scenes  are  stored  on  tape  in  an  8-bit  for¬ 
mat  (1  pixel  =  1  byte).  Because  of  this,  the  data  were 
scaled  to  range  from  0  to  255.  To  recover  radiance, 
a  multiplicative  factor  and  an  additive  factor  (AF)  were 
applied  to  the  data.  Each  scene  has  a  unique  MF  and 
AF.  These  factors  are  listed  in  Appendix  B.  Thirty 
scenes  were  available  for  analysis;  the  four  primary 
AFGL  scenes  used  in  this  analysis  were  101 R,  104R, 
105R,  and  I09R. 

The  third  set  of  data  used  in  this  study  was  from 
ERIM.  The  data  collection  instrument  is  a  line  scan¬ 


ner  flown  at  low  altitude.  Ten  spectral  bands  were 
available  with  the  ERIM  scanner.  At  most,  four  bands 
were  scanned  for  a  scene.  The  size  of  the  scenes  was 
variable.  The  FOV  of  the  scanner  also  varied  slight¬ 
ly.  Appendix  C  has  more  detailed  information  on  the 
ERIM  data. 

Six  ERIM  scenes  are  used  in  this  study.  These  are 
denoted  SOI,  S02,  S03,  S23,  S24,  and  F08.  There  were 
three  or  four  channels  for  each  of  these  scenes.  The 
first  three  of  these  scenes  are  from  Camp  A.P.  Hill, 
Virginia.  These  images  provide  morning,  afternoon, 
and  evening  data  of  one  scene.  S23  and  S24  are  of  Port 
Hueneme,  California.  S23  is  of  the  port  and  water; 
sunglint  is  added  to  this  in  S24.  The  last  image,  F08, 
is  of  the  mountains  by  Nellis  AFB. 

Table  3  compares  the  characteristics  of  the  BMAP, 
AFGL,  and  ERIM  data. 

Figure  1  shows  an  image  from  each  of  the  sets  of 
data  used  in  this  study.  Part  A  shows  middle- 
wavelength  BMAP  data  from  Bedford  (31  frames  of 
data  are  shown).  Part  B  of  the  figure  shows  the  AFGL 
scene  109R.  The  remainder  of  the  images  in  this  fig¬ 
ure  are  ERIM  images.  Part  C  is  S01CH2,  part  D  is 
S02CH2,  and  part  E  is  S03CH2.  These  three  are 
middle-wavelength  images  from  Camp  A.P.  Hill  at 
different  times  of  day.  Part  F  of  the  figure  is  S23CH3 
from  Port  Hueneme.  Finally,  part  G  of  the  figure  is 

F8CH 1 ,  a  mountain  scene  from  Nellis  AFB.  These  im¬ 
ages  will  be  shown  after  each  of  the  processing  steps 
to  show  the  effect  of  the  processing. 

Table  4  summarizes  the  data  used  in  this  study.  All 
of  the  images  used,  IFOV,  size,  spectral  band,  and 
descriptive  comments  are  included  in  this  table. 
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Table  3  Characteristics  of  instruments  used  to  coliect  background  data  used  in  this  study. 


Parameter 

BMAP  sensor 

AFGL  sensor 

ERIM  sensor 

Platform 

Ground 

NKC-135  @  35,000  ft 

C-47 

Viewing 

5“  elevation 

Typically  small 

90“  depression  (nadir)  ±45° 

geometry 

depression  angles 

scan  or  35°  depression  with 
±45°  scan 

Sensor 

3.9-4.8 

7.84-13.16 

various  bands: 

band 

7.6-11.3 

2.0-2.6 

(urn) 

3. 0-4.6 

Instantaneous 
field  of  view 
(mrad) 

Sensitivity 
NEI  (W/cm^) 


NEAT  CO 


Scene  size 


0.33  X  0.33 
(azimuth 
resolution 
=  0.096) 


Long  1.1  X  10  • 
Mid  6.1  X  10“'^ 


0.36  (H)  X  0.43  (V) 
(digitized  sample  size 
=  0.5  X  0.5) 


5.5  X  10  '^ 


Long  0.02 
Mid  0.01 


8  scans  mid-IR  and 
8  scans  long-1  R 
per  frame 

or 

16  scans  mid-IR  per 
frame  and  15  scans 
long-IR  per  frame 

373  pixels  per  scan 
line  oversampled  at 
3.3:1 


256  (H)  X  192 (V) 
pixels 


3.5- 3.9 
3.9-4.7 

4.5- 5. 5 
8.0-14.0 
9.0-11.4 

2.0  or  2.5  (circular) 


2.0-2.6  iim  4.91  x  10" 

3. 5- 3.9  tim  8.83  x  10' 
3.9-4.7  ftm  3.40  x  10" 

4.5- 5.5  ^xm  4.47  x  10' 

8.0-14.0  2.45  x  10' 

9.0-11.4  nm  4.63  x  10 

3. 5- 3.9  ,xm  3.0 
3.9-4.7  ^lm  0.25 

4.5- 5. 5  fim  0.2 
8.0-14.0  urn  0.1 
9.0-11.4  fxm  0.25 

855  or  646  pi.xels  (H); 
number  of  lines  per  scene 
varies  from  700  to  2000 
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Figure  1  Images  from  each  of  the  sets  of  data  used 
in  this  study,  (a)  BMAP  middle-wavelength  image  from 
Bedford,  (b)  AFGL  scene  I09  of  Monument  Valley,  (c) 
Morning  (scene  S01CH2).  (d)  Afternoon  (scene 
S02CH2).  (e)  Evening  (scene  S03CH2).  (f)  ERIM  image 
of  Port  Hueneme  (4.5-5.5  spectral  band)  (scene 
S23CH2).  (g)  ERIM  Image  of  mountains  by  Nellis  AFB 
(3.9-4.7  ^m  spectral  band)  (scene  F8CH1). 
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Table  4  Summary  of  data  used  in  this  study. 


Data 

Scan  IFOV  (H) 

Image  size 

source 

(mrad) 

(H) 

(V) 

Comments 

BMAP 

0.33 

373 

8 

Montauk  (50  frames); 

long-  and  mid-IR;  overcast  sky; 

sample  spacing:  0.1  mrad 

BMAP 

0.33 

372 

15 

Bedford  1  (61  frames) 

Long-IR;  hole  in  cloud  deck 

16 

Mid-IR;  hole  in  cloud  deck 

Sample  spacing:  0.1  mrad 

File  Description 


AFGL 

0.5 

256 

192 

lOlR 

clouds,  down  viewing 

104R 

terrain,  down  viewing 

105R 

horizon,  clouds,  terrain 

109R 

Monument  Valley,  clouds 

File 

Band 

Camp  A.P.  Hill 

ERIMI 

2.0 

855 

700 

SOICHI 

2.0-2.6 

Morning 

700 

S01CH2 

4. 5-5. 5 

m 

700 

S01CH3 

8.0-14.0 

M 

704 

S02CH1 

2.0-2.6 

Afternoon 

704 

S02CH2 

4.5-5.5 

704 

S02CH3 

8.0-14.0 

H 

700 

S03CH1 

Bad  data 

Evening 

700 

S03CH2 

4.5-5. 5 

M 

700 

S03CH3 

8.0-14.0 

File 

Band 

Port  Hueneme 

ER1M3 

2.5 

646 

1982 

S23CH1 

2.0-2.6 

Port/water 

1982 

S23CH2 

Bad  data 

M 

1982 

S23CH3 

4.5-5. 5 

1982 

S23CH4 

9.0-11.4 

1761 

S24CH 1 

2.0-2.6 

Port/water/sun 

1761 

S24CH2 

3.0-4. 2 

1761 

S24CH3 

4. 5-5. 5 

1761 

S24CH4 

9.0-11.4 

File 

Band 

Nellis  AFB 

ER1M5 

2.5 

646 

1541 

F08CH1 

3.9-4. 7 

Mountains 

1541 

I08CH2 

3. 5-3. 9 
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3.0  CHARACTERIZATION  OF  DATA 


3.1  Data  Artifacts  and  Corrections 

There  are  a  number  of  data  artifacts  and  charac¬ 
terization  uncertainties  that  limit  the  usefulness  of 
scene  data  to  background  investigations.  All  data  used 
in  this  investigation  have  one  or  more  of  the  follow¬ 
ing  unwanted  artifacts:  (a)  significant  sensor  noise,  (b) 
data  errors,  (c)  calibration  uncertainty,  (d)  imprecisely 
known  sensor  characteristics,  and  (e)  sensor-induced 
scene  artifacts. 

When  the  sensor  noise  masks  part  of  the  scene-clut¬ 
ter  power  spectrum,  the  high-spatial-frequency  com¬ 
ponent  of  the  background  clutter  is  lost  and  the  full 
spatial  resolution  of  the  collection  sensor  cannot  be 
used  in  evaluating  the  effect  of  clutter.  Data  errors 
can  include  poorly  performing  detectors,  recording  er¬ 
rors,  and  externally  produced  noise  (e.g.,  from  vibra¬ 
tion).  Such  errors  can  produce  target-like  signals  that 
significantly  degrade  usefulness  in  determining  the  ef¬ 
fectiveness  of  an  algorithm  that  uses  spatial  extent  as 
the  discriminant. 

Calibration  errors  include  radiometric  inaccuracy 
that  causes  a  significant  variation  in  the  predicted  per¬ 
formance  in  clutter  and  calibration  imbalance  between 
detectors  that  degrades  the  scene  in  the  cross-channel 
direction.  This  latter  effort  is  not  especially  signifi¬ 
cant  if  all  processing  is  done  using  data  from  the  col¬ 
lection  instrument  scan  direction  only.  Uncertainty  in 
the  clutter  processing  results  is  also  introduced  by  un¬ 
known  or  imprecise  description  of  the  instrument  and 
its  functioning.  Accurate  knowledge  of  the  spectral 
response,  and  optical  and  electrical  transfer  functions 
are  frequently  imprecise  or  unknown. 

Finally,  the  data  collection  geometry  can  adversely 
affect  the  data  by  introducing  aperture  vignetting  (i.e., 
lower  scene  variations)  or  geometric  distortion  (such 
that  the  pixel  footprint  on  the  scene  varies  significantly 
across  a  scan  line).  These  effects  are  particularly 
noteworthy  in  wide  FOV  instruments.  Another  geo¬ 
metric  effect  is  imprecise  or  changing  registration  be¬ 
tween  detectors  making  simultaneous  measurements 
in  multiple  spectral  bands.  Registration  between  spec¬ 
tral  bands  becomes  particularly  significant  if  multi¬ 
ple  colors  are  part  of  the  discrimination  processes. 

i  All  data  sets  used  in  this  investigation  had  artifacts 

and,  to  some  extent,  imprecise  knowledge  of  the  sen¬ 
sor  characteristics.  A  description  of  the  major  artifacts 
in  each  data  set  is  given  below  . 

‘  The  BMAP  data  had  some  defective  data  points 

caused  by  the  data  acquisition  system.  Software  creat¬ 


ed  at  the  Naval  Research  Laboratory  (NRL)  located 
the  bad  data  points  in  the  BMAP  images  and  correct¬ 
ed  them  by  linearly  interpolating  from  adjacent  good 
data  points.  These  interpolated  points  were  iden¬ 
tified  in  the  data  supplied  to  APL  by  assigning  them 
a  negative  value.  When  the  data  are  used,  these  nega¬ 
tive  values  are  changed  back  to  positive  values.  (For 
this  study,  we  did  not  need  to  keep  track  of  the  bad 
data  points.)  Correcting  these  data  recording  errors 
was  all  that  had  to  be  done  to  the  BMAP  data. 

When  the  BMAP  Montauk  images  were  first  re¬ 
ceived  at  APL  some  of  the  recording  errors  had  not 
been  corrected  by  the  NRL  software.  The  errors  were 
detected  after  the  images  were  processed  through  the 
clutter-suppression  algorithm.  (TTie  clutter-suppression 
algorithm  passed  many  of  the  bad  data  points  as  pos¬ 
sible  targets.)  NRL  then  sent  a  corrected  set  of  the  im¬ 
ages  to  APL.  The  performance  of  the  clutter-sup¬ 
pression  algorithm  was  much  better  on  these  updated 
images  than  on  the  original  images. 

Both  the  AFGL  and  some  ERIM  images  were  sig¬ 
nificantly  affected  by  noise  within  the  sensor.  This 
noise  masks  the  highest  frequencies  of  the  background 
clutter.  The  effect  of  this  noise  will  be  discussed  when 
the  power  spectral  densities  (PSD)  are  presented. 
Methods  of  reducing  the  effect  of  this  noise  will  also 
be  discussed  later  in  this  report. 

AFGL  images  were  formed  by  a  TV-like  raster  scan, 
with  alternate  lines  produced  by  separate  detector  ar¬ 
rays.  The  two  detector  arrays  produced  alternating 
lighter  and  darker  scan  lines  when  the  response  of  the 
two  arrays  had  drifted  out  of  calibration.  Reference 
14  describes  how  the  images  were  corrected.  First,  the 
mean  of  each  line  was  computed.  Half  of  the  differ¬ 
ence  of  the  mean  values  was  then  added  to  pixels  in 
the  darker  lines  and  half  the  difference  was  subtract¬ 
ed  from  pixels  in  the  lighter  lines.  All  of  the  AFGL 
scenes  shown  in  this  study  have  been  corrected. 

There  were  multiple  problems  with  the  ERIM  data. 
The  main  problem  was  that  images  obtained  with  a 
depression  angle  less  than  90°  suffered  from  vignet- 


'  K  A.  Sicinberg.  tiintinuliiin  ol  Sensor  Arii/acis  from  Infrared 
Dula.  NRI  Report  SHbl  (l)iv  II,  1984). 

'R  A  Steinberg  iind  M.  .1.  Mettugh.  .In  hrror  Detection  and 
Sinoothinti  Mitorithm  tor  Infrared  Data.  IHL'  API  TCi  l.tS.S 
(Apr  1986) 

"a.  >•  SaLreek.  \l(il  Infrared  lina«e  Restoration.  IHl'  API 
I  ll)(l|8\  I  (P9 
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ting  by  part  of  the  aircraft  fuselage.  Reference  15 
describes  how  this  problem  was  corrected.  A  down- 
scan  average  of  all  of  the  lines  in  the  scene  was  com¬ 
puted.  The  mean  and  standard  deviation  of  the  left 
and  right  edges  of  the  scene  were  then  adjusted  to  cre¬ 
ate  a  constant  mean  and  standard  deviation  for  the 
downscan  average.  This  restoration  improved  the  im¬ 
ages  but  some  bright  and  dark  bands  still  remained 
on  the  edges  of  the  restored  images. 

In  addition  to  the  vignetting  problem,  there  were 
some  other  problems  with  particular  scenes.  These 
problems  were  found  while  analyzing  the  scenes  dur¬ 
ing  the  study.  All  three  channels  of  SOI  had  10  bad 
scan  lines.  These  lines  were  found  because  they  had 
pixel  values  of  0  or  255.  No  other  lines  in  the  scene 
had  these  pixel  values.  This  artifact  (see  Fig.  Ic)  was 
corrected  by  simply  skipping  the  10  bad  lines. 

Channel  1  of  S03  was  not  processed.  ERIM  reports 
analyzing  these  data  never  included  this  scene  in  the 
analysis.  A  PSD  plot  of  the  scene  differed  from  PSDs 
of  all  of  the  other  scenes.  Most  likely,  the  sensor  was 
not  able  to  acquire  background  data  because  the 
signal-to-noise  ratio  was  too  low.  The  image  is  essen¬ 
tially  featureless.  This  scene  was  therefore  excluded 
from  the  study. 

Channel  2  of  S23  was  also  excluded  from  the  study. 
While  analyzing  this  image  (shown  in  Fig.  2)  some  un¬ 
expected  results  were  obtained.  Examination  of  the 
digital  data  of  this  scene  showed  that  the  pixel  values 
were  mostly  0  and  255.  Since  the  data  showed  that 
there  was  a  problem  with  this  image,  it  was  excluded. 

3.2  Radiance  Statistics 

Before  any  of  the  IR  data  were  processed  certain 
characteristics  of  the  data  were  needed.  A  histogram 
of  each  image  was  computed.  From  this  histogram  the 
maximum  and  minimum  in  the  scene  were  obtained. 
The  mean  and  variance  were  also  calculated  for  each 
image.  Table  5  lists  the  radiance  values  (mean,  vari¬ 
ance,  maximum,  and  minimum)  for  most  of  the  im¬ 
ages  used  in  this  investigation. 

The  radiance  statistics  are  needed  to  assess  the  per¬ 
formance  of  the  processing  techniques  used.  The  max¬ 
imum  and  minimum  irradiance  (with  mean  removed) 
will  be  compared  to  the  maximum  and  minimum  af¬ 
ter  processing.  This  comparison  provides  a  measure 
of  the  performance  of  the  processing  techniques. 


'‘a.  N.  Vavreck.  ERtM  tnfrarrd  Image  Restoration,  JHU/APL 
FID(I)«A.U4X)3. 


Figure  2  Image  showing  bad  data  in  ERIM  scene 
S23CH2. 


3.3  Power  Spectral  Density 

In  addition  to  calculating  the  moments  of  the  im¬ 
ages,  the  PSDs  were  also  computed  because  they  were 
needed  in  selecting  the  appropriate  Filters  for  the  im¬ 
ages.  Figure  3  shows  the  PSDs  of  the  BMAP  Mon- 
tauk  long-  and  mid-lR  bands.  Figure  4  shows  the  PSD 
of  one  of  the  AFGL  scenes.  Figure  5  shows  the  PSD 
of  an  ERIM  scene.  Appendix  D  has  the  PSDs  of  many 
of  the  remaining  scenes.  From  Figs.  3  through  5,  it 
is  seen  that  all  of  the  PSDs  have  approximately  the 
same  slope.  The  power  falls  off  as  the  inverse  of  the 
square  of  the  spatial  frequency  Such  a  fre¬ 

quency  dependence  is  often  observed  in  background 
data.  Figures  3  through  5  also  show  that  at  the  higher 
frequencies  all  of  the  PSDs  begin  to  flatten. 

The  level  at  which  a  PSD  flattens  is  the  noise  level 
of  the  system.  All  three  figures  show  the  flattening. 
These  Figures  and  the  additional  power  spectral  den¬ 
sity  data  given  in  Appendix  D  show  that  all  of  the 
AFGL  and  some  of  the  ERIM  PSDs  have  this  flat¬ 
tening.  BMAP  data  also  show  some  change  in  slope 
at  high  frequencies. 

The  flattening  of  BMAP  PSDs  occurs  because  of 
the  oversampling  of  the  scene.  The  BMAP  detectoi 
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Table  5  Statistics  of  background  data  in  units  of  W/cm^/sr. 


Image  Statistics 


Scene 

Spectral 

band 

(Min) 

Mean 

Variance 

Max 

Min 

A 

(max  -  min) 

BMAP-Montauk 
(long  IR) 

7.6-11.3 

1.93 

X 

10“^ 

1.43  X  10“’ 

2.04 

X 

10  ’ 

1.84  X 

10“’ 

0.20  X  10“’ 

BMAP-Montauk 
(mid  IR) 

3. 9-4.8 

7.35 

X 

10"^ 

2.17  X  10"'^ 

9.10 

X 

10“^ 

6.41  X 

10“’ 

2.69  X  10“’ 

BMAP-Bedford 
(long  IR) 

7.6-11.3 

1.66 

X 

10“^ 

5.09  X  10 

1.76 

X 

10“^ 

1.30  X 

10“’ 

0.46  X  10“’ 

BMAP-Bedford 
(mid  IR) 

3.9-4.8 

1.05 

X 

10'^ 

3.12  X  10“^ 

1.92 

X 

10“'' 

7.26  X 

10“’ 

1.19  X  10““ 

AFGL 

101 R 

7.84-13.16 

2.87 

X 

10  ^ 

5.01  X  10“* 

3.87 

X 

10“’ 

2.42  X 

10“’ 

1.45  X  10“’ 

AFGL 

I04R 

7.84-13.16 

3.81 

X 

10  ’ 

1.36  X  10  ’ 

4.43 

X 

10“’ 

2.74  X 

10“’ 

1.69  X  10  ’ 

AFGL 

I05R 

7.84-13.16 

3.01 

X 

10  ^ 

2.69  -  10  * 

3.42 

X 

10“’ 

2.44  X 

10“’ 

0.98  X  10“’ 

AFGL 

I09R 

7.84-13.16 

2.80 

X 

10'^ 

3.57  X  10  * 

3.46 

X 

10“* 

2.32  X 

10“’ 

1.14  X  10“’ 

ERIM 

F8CH1 

3. 9-4. 7 

5.2 

X 

10  ’ 

7.8  X  10“" 

7.2 

X 

10  ^ 

2.5  X 

10“’ 

4.7  X  10“’ 

ERIM 

F8CH2 

3.5-3.9 

1.2 

X 

10-’ 

2.0  X  10  “ " 

2.8 

X 

10  ’ 

2.3  X 

10“* 

2.8  X  10“’ 

ERIM 

F8CH3 

9.0-11.4 

1.8 

X 

10  ^ 

2.8  X  10“* 

2.3 

X 

10  ’ 

1.4  X 

10  ’ 

0.9  X  10  ’ 

ERIM 

SOICHI 

2.0-2.6 

3.6 

X 

10 

2.6  X  10  "* 

1.0 

X 

10  ^ 

0.0 

1.0  X  10  “ 

ERIM 
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Figure  3  Power  spectral  density  of  6MAP  Montauk 
Point  middle-wavelength  and  long-wavelength  back¬ 
ground  data. 


Figure  5  Power  spectral  density  of  ERIM  Port  Hue- 
neme  scene  in  2.0-2.6  /xm  spectral  band  (scene 
S23CH1). 


Figure  4  Power  spectral  density  of  AFGL  cloud  and 
terrain  background  scene  105. 


has  a  resolution  of  0.33  mrad,  giving  a  maximum  spa¬ 
tial  frequency  of  about  1.5  cycles/mrad.  However,  the 
scene  is  sampled  at  0.1  mrad  resolution,  giving  a  max¬ 


imum  Nyquist  frequency  of  5  cycles/mrad.  Since  the 
detector  limits  the  frequency  to  1 .5  cycles/mrad,  data 
between  1.5  cycles/mrad  and  5  cycles/mrad  are  prin¬ 
cipally  sensor  noise. 

All  of  the  AFGL  PSDs  show  the  flattening  because 
at  high  spatial  frequencies  the  power  of  the  noise  with¬ 
in  the  sensor  is  greater  than  the  power  from  the  scene. 
Without  the  sensor  noise,  the  PSDs  are  expected  to 
maintain  the  \/k^  slope  up  to  1 .0  cycle/mrad.  Sever¬ 
al  AFGL  scenes  with  higher  spatial  resolution  (-0.1 
mrad)  were  also  analyzed  and  show  that  the  \/k‘ 
slope  continues  to  at  least  this  spatial  frequency  (see 
Appendix  D). 

Table  6  lists  the  levels  at  which  the  four  PSDs  flat¬ 
ten  and  the  frequencies  where  they  begin  to  flatten. 
These  values  are  also  listed  for  some  of  the  PSDs  giv¬ 
en  in  Appendix  D.  Also  listed  in  Table  6  are  the  cut¬ 
off  frequencies  of  the  sensors.  The  roll-off  of  the 
PSDs,  the  noise  level,  and  the  cutoff  frequency  are 
used  in  the  analysis. 

The  point-source  characteristic  of  the  BMAP  data 
is  included  in  Appendix  A.  Since  the  BMAP  scenes 
are  oversampled,  a  point  target  in  the  scene  will  spread 
into  many  pixels  in  the  resulting  image.  The  Fourier 
transform  of  the  point-source  characteristic  was  com¬ 
puted  for  both  the  long-  and  mid-IR  bands.  Figure 
6  shows  the  Fourier  transform  of  the  point-source 
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Table  6  PSD  noise  level. 


Power  spectral  density  System  cutoff 

at  noise  level  Flattening  frequency  frequency 

Scene  ((W/cm’/sr)^/(cycles/mrad))  (cycles/mrad)  (cycles/mrad) 


BMAP  (long  IR) 

1.5 

X 

10 

\y 

3.0 

1.5 

BMAP  (mid  IR) 

3.0 

X 

10 

15 

3.0 

1.5 

AFGL  (101 R) 

3.3 

X 

10 

10 

0.35 

1.0 

AFGL  (I04R) 

2.0 

X 

10 

9 

0.35 

1.0 

AFGL  (I05R1 

3.0 

X 

10 

10 

0.30 

1.0 

AFGL  (I09R) 

4.0 

X 

10 

10 

0.30 

1.0 

ERIM  (SOICHI) 

6.0 

X 

10 

i: 

0.15 

0.25 

ERIM  (S01CH2) 

1.5 

X 

10 

12 

0.17 

0.25 

ERIM  (S23CH1) 

2.0 

X 

10 

it 

0.15 

0.20 

Figure  6  Fourier  transform  of  point  source  as  seen 
by  the  BMAP  instrument. 


characteristic  of  the  long-IR  band.  The  plot  looks  like 
a  low-pass  filter  with  a  cutoff  at  about  1.5  cy¬ 
cles/mrad. 

Point-source  characteristics  of  the  AFGL  data  were 
not  available.  Some  ERIM  frequency  response  data 
are  included  in  Appendix  C.  Since  there  was  no  over- 
sampling  in  either  of  these  data  sets,  it  was  assumed 
that  the  point-source  target  is  an  impulse.  The  PSD 
of  an  impulse  is  a  flat  spectrum  with  a  cutoff  at  the 
Nyquist  frequency.  Therefore,  while  the  PSD  of  the 
BMAP  point-source  characteristic  looks  like  a  real  low- 
pass  filter,  the  PSDs  of  the  AFGL  and  ERIM  point- 
source  characteristics  are  assumed  to  be  ideal  low-pass 
filters. 
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4.0  SCENE  CLUTTER  PROCE.SSING 


4.1  Approach 

The  first  technique  examined  uses  the  PSD  of  the 
scene  to  give  an  estimate  of  the  required  target  signal 
for  detection.  A  signal-to-noise  criterion  was  selected 
to  produce  a  nominal  false-alarm  level.  The  target  is 
modeled  as  a  point  source  with  unknown  irradiance 
in  the  spectral  band.  A  matched  filler  is  then  applied 
to  the  signal  and  noise  spectra  to  estimate  the  required 
target  irradiance  for  detection.  The  target  irradiance 
required  for  detection  and  the  false-alarm  rate  are  the 
primary  measures  of  performance,  both  in  this  study 
and  for  tactical  sensors.  This  technique  is  well- 
established  and  often  used  to  estimate  IR  sensor  per¬ 
formance. 

A  matched-filter  approach  offers  a  number  of  ad¬ 
vantages.  Foremost  is  that  this  approach  allows  one 
to  make  performance  estimates  for  a  scanning  sensor 
in  the  presence  of  background  clutter  rather  than  on 
sensor  noise  alone.  Since  the  matched-filter  technique 
gives  the  maximum  achievable  performance,  this  meth¬ 
od  sets  a  limit  (or  standard)  on  sensor  detection  capa¬ 
bility.  The  matched-filter  approach  can  be  used  to 
compare  clutter-limited  performance  in  different 
scenes  and  in  different  spectral  bands  when  the  same 
scene  is  measured  in  multiple  colors. 

The  matched-filter  approach  to  predicting  clutter- 
limited  performance  has  well-known  disadvantages. 
Results  of  the  method  are  incorrect  because  real  scenes 
do  not  meet  the  assumptions  of  the  matched-filter 
analysis  (see  below).  Even  so,  a  processing  filter  that 
is  approximately  "matched”  to  the  signal  and  ciutter 
spectrum  approximately  minimizes  the  clutter  power 
into  the  sensor.  Hence  the  matched-filler  method  "de¬ 
signs”  a  real  filter  that  can  be  applied  directly  to  the 
raw  scene  data.  Directly  filtering  the  raw  scene  pro¬ 
vides  a  measure  of  false-alarm  rate  as  a  function  of 
detection  threshold.  The  threshold  can  then  be  used 
to  estimate  the  required  target  irradiance. 

Therefore,  the  second  approach  is  directly  filtering 
the  raw  scene  with  a  digital  filter  that  approximates 
a  matched  filter.  Raw  data  is  processed  with  a  first- 
order  high-pass  filter  that  approximately  "whitens" 
the  clutter  spectrum.  Results  of  this  filtering  process 
can  be  compared  to  that  of  the  matched-filter  analy¬ 
sis.  Several  methods  of  comparison  are  possible:  one 

'H  SamucKcin.  "Infrarccl  SiMcms;  I  I  xprc'Muii'  tor  Sil’ii.iI 
.ind  Background  Induced  Noise  vsith  Space  I  illers."  //  /  /  frun\ 
tvnop.  Hciiroit  Swt  4KS-7.  2'  lld"l| 

R  Yoshilani,  "Signal  Processing  Ilieors  lor  I  inearlv  Scanning 
Rectangular  IJeleclor."  f’r,,,  IKIS  18.  1 1'  UK  tIV'Tl 


is  to  determine  the  threshold  required  to  give  no  false 
alarms  in  the  filtered  data  and  then  compare  the  results 
with  the  threshold  that  would  be  required  to  give  no 
false  alarms  (e.g.,  false  alarm  probability  of  <0.5) 
in  the  matched-filter  analysis. 

Thresholding  the  filtered  data  at  various  levels  gives 
the  number  of  false  alarms  (i.e.,  false-alarm  rate)  as 
a  function  of  threshold  (hence,  target  signal).  After 
filtering,  the  scene  has  the  mean  background  level  re¬ 
moved  and  large,  low-frequency  features  are  sup¬ 
pressed.  Unfortunately,  large  features  with  sharp  edges 
will  also  produce  a  significant  filter  output  that  is  in¬ 
distinguishable  from  that  produced  by  a  target.  Real 
scenes  have  a  number  of  these  boundaries  such  as 
cloud  edges  against  sky,  terrain,  or  other  cloud  lay¬ 
ers,  terrain  crests  and  valleys,  and  boundaries  between 
woods  and  fields. 

Further  discrimination  is  obtained  by  examining  the 
differences  in  high-pass  filter  response  to  point  ob¬ 
jects  and  to  extended  objects.  Since  the  high-pass  fil¬ 
ter  acts  as  a  differentiator,  a  filtered  positive  point 
feature  produces  a  positive  pulse  immediately  followed 
by  an  approximately  equal-amplitude  negative  pulse. 
(The  reverse  is  true  for  a  negative  point  feature.)  On 
the  other  hand,  a  sharp  edge  will  only  produce  a  posi¬ 
tive  (or  negative)  pulse,  not  a  pair  of  pulses.  There¬ 
fore  a  pulse-width  discriminator  can  be  developed  to 
separate  small  point  features  from  edges. 

Implementation  of  the  pulse-width  discriminator  re¬ 
quires  that  a  negative  crossing  of  a  specially  set  con¬ 
firmation  threshold  will  follow  a  positive  crossing  of 
the  detection  threshold.  Similarly,  crossings  of  the 
negative  detection  threshold  (i.e.,  to  detect  a  negative 
contrast  target)  must  be  followed  by  a  crossing  of  the 
positive  confirmation  threshold.  The  confirmation 
thresholds  are  only  active  for  a  short  time  after  a  cross¬ 
ing  of  the  detection  threshold. 

This  pulse-width  discriminator  rejects  sharp  edges 
and  therefore  enhances  the  discrimination  provided 
by  the  high-pass  digital  filter.  Performance  is  again 
measured  in  terms  of  number  of  detections  as  a  func¬ 
tion  detection  threshold  leve'  (i.e.,  false-alarm  rate). 
However,  the  high-pass  filter  plus  pulse-width  discrimi¬ 
nator  is  still  not  sufficient  to  eliminate  extended  scene 
features  totally.  I  inear  features  such  as  roads  and  sun- 
enhanced  cloud  edges  will  pass  the  processing  tests. 
If  the  sensor  has  multiple  parallel  detection  channels, 
near  simultaneous  detections  in  adjacent  channels  can 
be  used  to  reject  features  with  extent  in  the  cross-scan 
direction. 
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This  post-detection  adjacent-channel  logic  will  al¬ 
low  discrimination  against  linear  features  whose  length 
extends  across  three  or  four  contiguous  channels.  The 
combination  of  high-pass  filtering,  pulse-width  dis¬ 
crimination,  and  adjacent-channel  logic  provides  a 
complete  set  of  spatial-discrimination  techniques  that 
allow  separation  of  point  features  from  larger  features 
with  significant  in-scan  or  cross-scan  extent  (or  both). 
The  nature  of  the  scanning  geometry  provides  better 
discrimination  in  the  in-scan  direction  because  of  over- 
sampling. 

Even  implementation  of  the  above  combination  of 
techniques  cannot  eliminate  all  potential  false  alarms. 
Background  features  whose  spatial  extent  is  <2  de¬ 
tector  widths  (or  pixels  or  resolution  cells)  in  the  scan 
direction  and  <3  detector  widths  in  the  cross-scan 
direction  cannot  be  distinguished  from  targets.  As  sen¬ 
sor  sensitivity  is  improved,  the  number  of  these  small 
features  that  pass  all  point  source  discrimination  tests 
will  increase  and  clutter  will  limit  the  sensor  perfor¬ 
mance.  When  this  fundamental  limit  is  encountered, 
other  discrimination  techniques  such  as  improved  spa¬ 
tial  resolution  (unlikely  in  a  missile  seeker)  or  color 
(multiple-spectral-band)  tests  must  be  investigated. 
Such  an  investigation  requires  a  database  beyond  the 
scope  of  this  investigation. 

The  following  sections  describe  the  analyses  and 
results  of  these  discrimination  techniques  as  applied 
to  the  available  background  scenes. 


4.2  Matched-Filter  Performance  Analysis 

Using  the  PSDs  of  the  .AFGL,  BMAP,  and  ERIM 
images,  a  matched  filter  was  applied  to  the  data. 

The  signal  was  assumed  to  be  an  impulse  (flat  spec¬ 
trum  to  the  cutoff  frequency),  and  the  cutoff  frequen- 
c>  <i.e.,  spatial  resolution)  was  varied  up  to  the  sen.sor 
resolution-limited  Nyquist  maximum  value  to  deter¬ 
mine  the  effect  of  spatial  resolution.  A  minimum  tar¬ 
get  irradiance  required  for  detection  was  determined 
as  a  function  of  the  cutoff  frequency  assuming  a  tar¬ 
get  irradiance  to  rms  clutter  noise  irradiance  of  \50. 
Hence,  the  detection  threshold  is  7.07  times  the  rms 
clutter  noise.  Assuming  that  the  noise  is  Gaussian,  the 
false-alarm  probability  is  0.5  after  about  10'‘  in¬ 
dependent  looks  (resolution  cell). 


M  Schrocdcr.  .1.  H.  Schumnicrs.  B.  P.  Sandlord.  and  \V.  .1. 
rropl  .  "Intrarcd  t  loud  Backiirouiids."  in  Proc.  Trt-Svn  icc  hi- 
Irured  Hai  kgrotitui  Syiitpositim .  .M  til  -  rR-S4-0094  (Mar  16, 
19S4). 

'1  .  M.  Howscr  and  V\  .  .1.  Tropl,  ‘InrcMigalion  ol  Infrared 
Background  (  Inner."  in  Research  and  Pevelopinenl  Programs 
al  The  Johns  Hopkins  i  niversily  Applied  Physics  I  ahoralory. 
IHl  -XPl  yuarierly  Report  RQR  8.S-1  (.lan-Mar  19X5). 


From  the  above  discussion,  the  signal-to-noise  ra¬ 
tio  (SNR)  is  given  as 


SNR 


|s(0)|- 

0\ 


(1) 


where  s(0)  is  the  peak  signal  output  and  is  the 
variance  of  the  noise  output.  For  a  system  with  a  fil¬ 
ter  transfer  function  H{J)  the  peak  signal  and  noise 
variance  are 


5(0) 


H(J)  s(f,  df 


(2a) 


and 


N(f)  df 


(2b) 


where  s(f)  is  the  signal  spectrum  and  N(f)  is  the  noise 
(internal  plus  scene)  PSD.  For  a  matched  filter 


SNR 


'---i 


|5(/)|- 

N(f) 


df 


(3) 


and  when  the  target  is  an  impulse,  the  signal  spectrum 
(for  AFGL  and  ERIM  data)  is  assumed  to  be  gener¬ 
ated  from  an  impulse: 


S(k)  =  E,  Xx  (4) 

where  E,  is  the  target  irradiance  (in  W/cm")  and  Xc 
is  the  sample  interval.  Table  7  gives  seeker  perfor¬ 
mance  results  in  terms  of  required  target  irradiance 
for  detection  using  the  above  method  for  all  available 
AFGL  scenes.’*’  Three  cases  are  given.  The  first  uses 
the  PSDs  of  the  original  data  (background-plus-noise 
case  that  includes  the  collection  scenes).  The  second 
performance  result  is  the  required  target  irradiance  if 
only  sensor  noise  from  the  forward-looking  IR  (FLIR) 
were  present  (FLIR-noise-only  case).  The  AFGL  FLIR 
internal  noise  is  approximated  by  a  flat  PSD  with  a 
magnitude  equal  to  the  high  spatial  frequency  portion 
of  the  PSD.  Since  the  AFGL  data  always  show  sen¬ 
sor  noise  at  high  spatial  frequencies,  this  internal  noise 
estimate  gives  the  sensor  limitation.  Finally,  the  third 
column  (background-only  case)  gives  an  estimate  of 
the  background  clutter  limitation,  assuming  that  the 
clutter  PSD  continues  to  decrease  as  \/k'.  Averages 

"VV.  ,1.  Tropl,  A.  N.  Vavrcck,  B.  .1.  .Saiidford,  ,1.  H.  Scliunimcrs. 
and  .1.  Schrocdcr,  "Inlrcired  C  loud  Backgrounds  and  Sensor  Per- 
lorinancc,"  in  Presenlations  ai  the  Second  Tri-Service  Cloud 
Modeling  Workshop.  ID,\  Report  M-9,  I,  pp.  199-248  (.Aug  I, 
1984). 
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Table  7  Required  target  irradiance  for  matched  filter  detection  in  several  AFGL  cloud  scenes. 

Required  Target  Irradiance  (10  "  W/cm^) 


Background 

FLIR 

Background 

Image  IFOV 

Scene  description 

plus  noise 

noise  only 

only 

101  0.5  mrad 

Down  viewing  cloud 

3.98 

2.96 

1.81 

102 

Down  viewing  and  water 

3.46 

2.67 

1.53 

103 

Down  viewing  and  terrain 

8.36 

6.46 

4.80 

104 

Down  viewing  and  terrain 

8.67 

6.82 

4.67 

105 

Clouds  and  horizon 

3.53 

2.73 

1.73 

106 

Cloud  deck 

4.97 

3.68 

2.95 

107 

Cloud  structure 

5.53 

4.22 

3.70 

108 

Clouds  and  horizon 

9.81 

8.14 

2.94 

109 

Monument  Vailey 

4.17 

3.21 

2.34 

no 

Monument  Valley 

4.09 

3.12 

2.23 

Ill 

Close  clouds 

3.39 

2.71 

1.50 

112 

Close  clouds 

3.38 

2.72 

1.63 

113 

Clouds  and  horizon 

5.29 

4.36 

2.62 

114 

Clouds  and  horizon 

5.75 

4.65 

3.01 

II6 

Horizon  and  clouds 

6.20 

4.62 

3.26 

117 

Clouds  and  terrain 

5.20 

4.35 

2.07 

118 

Clouds  and  horizon 

5.95 

4.57 

3.02 

119 

Clouds  and  horizon 

5.64 

4.34 

3.36 

120 

Down  viewing  clouds 

3.32 

2.69 

1.40 

131 

Clouds  below  horizon 

4.54 

3.55 

3.23 

132 

Down  viewing  clouds 

4.96 

3.86 

3.70 

21 -scene,  0.5-mrad  average 

5.23 

4.12 

2.74 

121  0.1  mrad 

Structured  clouds 

0.147 

0.110 

0.070 

122 

Structured  clouds 

0.152 

0.107 

0.083 

123 

Structured  clouds 

0.175 

0.131 

0.106 

124 

Structured  clouds 

0.150 

0.122 

0.072 

125 

Structured  clouds 

0.171 

0.134 

0.084 

126 

Structured  clouds 

0.154 

0.118 

0.077 

127 

Structured  clouds 

0.165 

0.125 

0.084 

128 

Structured  clouds 

0.160 

0.127 

0.077 

129 

Structured  clouds 

0.144 

0.114 

0.063 

9-scene,  0.1 -mrad  average 

0.157 

0.121 

0.080 

of  the  results  for  21  AFGL  scenes  sampled  at  0.5  mrad 
and  for  9  AFGL  scenes  sampled  at  0. 1  mrad  are  given. 

Similar  matched-filter  estimates  were  made  using 
BMAP  Montauk  Point  and  Bedford  No.  1  data.  Since 
scanning  seeker  performance  is  dependent  on  resolu¬ 
tion  (i.e.,  IFOV),  a  matched-filter  analysis  with  a  vari¬ 
able  seeker  IFOV  was  conducted.  Results  for  AFGL 
and  BMAP  cloud  scenes  are  given  in  Fig.  7.  The  tar¬ 
get  irradiance  required  for  detection  is  proportional 


to  the  sensor  instantaneous  solid-angle  FOV  (square 
of  linear  FOV  dimension).  AFGL  results  for  the  21 
scenes  with  0.5  mrad  resolution  and  9  scenes  with  0. 1 
mrad  resolution  are  shown  as  a  shaded  range  of  values. 
Results  of  two  BMAP  scenes  are  shown  for  compari¬ 
son.  Long-wavelength  BMAP  data  are  shown  as  thick 
solid  lines  and  middle-wavelength  results  as  thick 
dashed  lines.  The  Montauk  Point  scene  is  overcast  and 
has  low  contrast.  The  Bedford  No.  1  is  a  hole  in  the 
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Figure  7  Required  target  irradiance  for  detection  as  a  function  of  seeker  resolution  for  AFGL  and  BMAP  data. 


clouds  near  the  sun  and  has  significantly  higher  con¬ 
trasts.  Note  that  long-wavelength  clutter-limited  results 
from  both  AFGL  and  BMAP  data  are  comparable, 
although  the  overcast  sky  at  Montauk  Point  gives  low¬ 
er  clutter  results. 

A  typical  middle-wavelength  seeker  will  have  a  spa¬ 
tial  resolution  of  O.S  to  1 .0  mrad  while  the  resolution 
of  a  long-wavelength  seeker  will  be  1  to  2  mrad.  Com¬ 
bining  the  difference  in  minimum  detectable  target  sig¬ 
nature  and  resolution  results  in  approximately  two 
orders  of  magnitude  difference  between  the  clutter 
noise  of  the  middle-  and  long-wavelength  bands. 

An  ideal  bandpass  filter  was  also  applied  to  the 
Montauk  data  using  the  PSDs.  ^  The  response  of  an 
ideal  bandpass  filter  is  zero  up  to  the  cut-on  of  the 
filter.  After  the  cut-on,  the  filter  response  is  unity  up 
to  the  cutoff.  The  cutoff  of  this  bandpass  was  chosen 
to  be  the  sensor-limited  maximum  frequency  (1 .5  cy- 
cles/mrad).  The  cut-on  of  the  filter  was  varied  from 
just  above  DC  to  the  sensor-limited  maximum  fre¬ 
quency.  The  minimum  detectable  signal  (or  required 
target  irradiance)  was  computed  as  a  function  of  cut- 


on  frequency,  assuming  a  signal-to-noise  ratio  of  7.07. 
Figure  8  shows  a  comparison  of  the  matched-filter  per¬ 
formance  versus  the  ideal  bandpass  filter  performance. 
Figure  8a  is  from  the  mid-IR  data  and  8b  is  from  the 
long-IR  data.  Both  parts  of  the  figure  show  that  the 
matched  filter  requires  a  lower  target  irradiance  than 
the  ideal  bandpass  filler.  However,  the  ideal  bandpass 
filter  performance  approaches  the  matched-filter  per¬ 
formance  for  a  bandpass  cut -on  that  is  approximate¬ 
ly  one-half  of  the  cutoff  frequency  of  the  system. 

Using  the  PSDs  of  some  of  the  AFGL  scenes,  an 
ideal  high-pass  filter  was  applied  to  AFGL  images. 

A  high-pass  filter  was  used  for  the  AFGL  data  because 
the  scene  was  sampled  at  the  Nyquist  rate.  Therefore, 
the  sensor-limited  maximum  frequency  equaled  the 
sampling  frequency,  and  the  higher  frequencies  did 
not  need  to  be  rejected.  The  results  of  the  procedure 
were  the  same.  The  matched  filter  performed  belter 
than  the  ideal  high-pass  filter,  but  the  performance 
of  the  high-pass  filter  again  approached  the  matched- 
filter  performance  for  a  cut-on  frequency  of  about 
one-half  the  cutoff  frequency  of  the  system. 
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Figure  8  Comparison  of  required  target  irradiance  for 
ideal  bandpass  fiiters  to  the  matched-filter  result  us¬ 
ing  BMAP  data,  (a)  Required  target  irradiance  for  BMAP 
mid-IR  images,  (b)  Required  target  irradiance  for  BMAP 
long-IR  images. 


The  matched  filter  and  ideal  high-pass  or  bandpass 
filter  results  give  a  first-order  approximation  to  seek¬ 
er  performance.  However,  these  results  assume  that 
the  scene  is  Gaussian  and  stationary.  The  stationary 
assumption  requires  that  the  power  spectrum  and 
statistics  of  the  scene  are  independent  of  the  location 
within  the  scene.  If  the  PSD  is  dependent  on  the  lo¬ 
cation  within  the  scene,  then  the  scene  average  PSD 
is  not  a  valid  representation.  Also,  the  PSD  charac¬ 
terization  of  the  scene  only  includes  a  measure  of  the 


scene  mean  and  variance,  not  higher  moments.  If  the 
scene  is  Gaussian  the  mean  and  variance  of  the  scene 
totally  describe  the  scene  statistics.  However,  if  the 
scene  is  not  Gaussian  (as  is  usually  the  case),  higher 
moments  are  needed. 

Since  real  scenes  violate  the  Gaussian  and  station¬ 
ary  assumptions,  actual  performance  (e.g.,  true-false 
alarm  rate)  must  be  determined  by  direct  processing 
of  the  raw  background  data.  The  following  sections 
of  the  report  describe  a  procedure  used  to  process  the 
raw  background  images. 

4.3  Digital-Filter  Performance  Analysis 

In  theory,  a  matched  filter  applied  to  the  IR  data 
would  maximize  the  signal-to-noise  ratio  of  a  point 
target  in  the  scene.  Since  the  assumptions  on  which 
matched-filter  theory  is  based  do  not  apply  to  IR  clut¬ 
ter,  direct  filtering  of  the  raw  scene  data  with  a  digi¬ 
tal  filter  that  approximates  a  matched  filter  can  be  used 
to  estimate  detection  system  performance. 

Since  the  signal  (point-source)  spectrum  is  approx¬ 
imately  white,  the  input  clutter  is  whitened  by  a  digi¬ 
tal  high-pass  filter  that  approximately  whitens  the 
clutter  and  therefore  represents  a  matched  filter.^' 

Since  the  BMAP  frequency  data  extend  beyond  the 
signal  (point-source)  spectrum,  an  additional  low-pass 
filter  is  needed  to  eliminate  sensor  noise.  The  addi¬ 
tional  filter  was  applied  to  the  BMAP  data  to  elimi¬ 
nate  data  lying  between  1 .5  and  5.0  cycles/mrad  spatial 
frequency. 

It  was  mentioned  earlier  that  the  PSDs  of  the  point- 
source  characteristics  of  the  AFGL  and  ERIM  data 
collection  instruments  were  close  to  flat  spectrums  be¬ 
cause  the  data  were  taken  at  the  Nyquist  frequency; 
therefore,  no  filtering  beyond  the  clutter-matched  filter 
was  needed  with  the  AFGL  and  ERIM  data. 

The  matched  filters  needed  for  these  images  were 
applied  to  the  PSDs  shown  in  Figs.  3  to  5.  A  Butter- 
worth  digital  filter  was  selected  to  evaluate  the  IR  im¬ 
ages.  “  A  first-order  high-pass  filter  was  used  because 
it  approximately  matched  the  rolloff  of  the  back¬ 
ground  PSDs  and  therefore  approximated  the  matched 
filter  of  the  previous  section.  PSDs  of  the  AFGL  and 
ERIM  filtered  images  are  approximately  flat.  The 
BMAP  PSDs  were  not  completely  flat  because  of  the 
second-order,  low-pass  filter  that  was  also  applied. 

The  cut-on  of  the  filters  was  chosen  on  the  basis 
of  the  point  where  the  PSDs  of  the  images  began  to 

‘'L.  M.  Howser,  Analysis  of  Background  Measurement  and 
Analysis  Program  (BMAP)  IR  Images,  JHU/APL 
FIF(2)85-U-035,  Rev  1. 

i’G.  L.  Turin,  "An  Introduction  to  Matched  Filters,"  IEEE 
Trans.  Information  Theory  IT  6,  311-329  (Jun  1960). 
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flatten.  Table  6  lists  the  system  cutoff  frequencies  for 
the  data  and  the  frequencies  at  which  the  PSDs  begin 
to  flatten  for  a  number  of  scenes.  These  frequencies 
were  used  to  choose  the  cut-on  of  the  highpass  filters. 
For  the  BMAP  data,  the  system  Nyquist  sampling  spa¬ 
tial  frequency  of  1.5  cycles/mrad  was  chosen  as  the 
cut-on  frequency.  Since  the  AFGL  image  flattened  at 
0.3  cycles/mrad,  0.25  was  chosen  as  the  filter  cut-on 
frequency.  (Since  the  specified  cut-on  frequency  is  3 
dB  down  for  a  first-order  high-pass  filter,  by  choos¬ 
ing  the  cut-on  frequency  to  be  slightly  less  than  the 


frequency  at  which  the  PSD  is  flat,  the  output  PSD 
should  be  approximately  flat.)  The  ERIM  cut-on  fre¬ 
quencies  were  chosen  for  the  same  reasons.  That  is, 
knowing  the  frequencies  where  the  input  PSDs  flat¬ 
tened,  the  cut-on  frequencies  were  chosen  to  whiten 
the  output  PSD. 

Table  8  lists  the  filters  and  the  corner  frequencies 
used.  Figures  9  through  1 1  show  the  filters  used  in  this 
analysis.  Figure  9  shows  the  BMAP  bandpass  filter. 
Fig.  10  shows  the  AFGL  high-pass  filter,  and  Fig.  1 1 
shows  the  ERIM  high-pass  filter. 


Table  8  Digital  filters  used  in  the  performance  analysis. 


Data 

source 

HFOV 

(mrad) 

Low-pass 

filter 

High-pass 

filter 

Corner 

(cycles/mrad) 

BMAP 

0.33 

2nd  order 

1st  order 

1.5 

AFGL 

0.5 

None 

1st  order 

0.25 

ERIM 

2.0 

None 

1st  order 

0.1 

ERIM 

2.5 

None 

1st  order 

0.07 

Figure  9  BMAP  bandpass  filter:  second-order  low- 
pass,  and  first-order  high-pass  (corner  =  1.5  cy¬ 
cles/mrad). 


Figure  10  AFGL  first-order  high-pass  filter  (cut-on  = 
0.25  cycle/mrad). 
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Figurell  ERIM  first-order  high-pass  filter  (cut-on  = 
0.1  cycle/mrad). 


Figure  12  Power  spectral  density  of  BMAP  Montauk 
Point  long-wavelength  background  data  after  bandpass 
filter. 


Figures  12  through  14  show  the  PSDs  of  three  im¬ 
ages  after  filtering.  Figure  12  shows  the  PSD  of  the 
BMAP  Montauk  long-IR  data  after  filtering  with  the 
combination  of  a  second-order,  low-pass  filter  and  a 
first-order,  high-pass  filter.  This  PSD  can  be  compared 
to  the  long-IR  PSD  in  Fig.  3  and  the  PSD  of  the  point- 
source  characteristic  in  Fig.  6.  The  comparison  with 
Fig.  6  shows  that  the  filtered  PSD  does  look  like  the 
Fourier  transform  of  the  point-source  characteristic, 
as  expected  after  applying  the  matched  filter.  Figures 
13  and  14  show  the  PSDs  of  the  AFGL  and  ERIM 
images.  These  PSDs  can  be  compared  to  the  PSDs  in 
Figs.  4  and  5.  The  two  figures  show  that  the  high-pass 
filters  did  approximately  flatten  the  background  PSDs. 
Appendix  E  shows  the  PSDs  of  the  remaining  images. 

The  reason  for  applying  the  filtering  was  to  max¬ 
imize  the  signal-to-noise  ratio  by  minimizing  the  clut¬ 
ter.  Flattening  the  PSD  of  the  background  image 
reduces  the  energy  of  the  low-frequency  objects  and 
does  not  affect  the  high-frequency  objects.  Point  tar¬ 
get  signals,  which  are  predominately  of  high-frequency 
content,  therefore  remain  almost  unaffected. 

Figures  15  and  16  demonstrate  the  effect  of  the 
filtering.  Figure  15a  shows  line  471  of  ERIM  scene 
SOICHI.  This  line  contains  no  point  targets.  Figure 
15b  shows  the  same  line  after  filtering,  now  all  of  the 
low-frequency  features  in  the  line  have  been  attenu¬ 
ated.  Figure  16a  shows  line  471  from  ERIM  scene 


Figure  13  Power  spectral  density  of  AFGL  cloud  and 
terrain  background  scene  I05  after  high-pass  filter. 


S01CH3.  This  line  has  some  possible  point  targets  pres¬ 
ent.  Figure  16b  shows  the  same  line  after  filtering; 
most  of  the  image  features  have  been  attenuated,  but 
the  possible  point  targets  are  almost  unaffected  by  the 
filter.  Figure  16c  shows  an  enlargement  of  the  area 
around  the  possible  point  targets,  giving  a  better  view 
of  the  filtered  output. 

The  previous  figures  have  shown  that  the  filtering 
applied  to  the  background  images  reduce  the  ampli¬ 
tude  of  image  features  when  possible  targets  are  not 
present.  Also,  the  mean  of  all  of  the  images  are  re¬ 
moved  when  high-pass  filters  are  applied. 
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Table  9  shows  the  maximum  and  minimum  radi¬ 
ance  of  the  scenes  used  in  this  study,  after  filtering. 
The  values  can  be  compared  to  the  values  in  Table  5. 
A  comparison  of  the  tables  shows  that  except  for  two 
scenes,  the  difference  between  the  maximum  and  mini¬ 
mum  scene  radiance  values  before  filtering  is  greater 
than  the  difference  after  filtering.  Figures  15  and  16 
show  that  this  result  is  expected. 


Figure  14  Power  spectral  density  of  ERIM  Port  Hue- 
neme  scene  in  2.0-2.6  spectral  band  (scene 
S23CH1)  after  high-pass  filter. 


0  300  600  900 

Pixel 


Figure  15  Line  471  of  ERIM  Camp  A.P.  Hill  morning 
scene  in  2.0-2.6  nm  spectral  band  (scene  S01CH1).  (a) 
Unprocessed  line  of  scene,  (b)  Line  of  scene  after  high- 
pass  filter. 


Figure  16  Line  471  of  ERIM  Camp  A.P.  Hill  morning 
scene  in  8.0-14.0  urn  spectral  band  (scene  S01CH3).  (a) 
Unprocessed  line  of  scene  with  data  artifact,  (b)  Line 
of  scene  after  high-pass  filter,  (c)  Expanded  line  of 
scene  after  high-pass  filter. 
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Table  9  Statistics  of  background  data  after  filtering  in  units  of  W/cm^/sr. 


Scene 

Maximum 

Minimum 

A 

(max  - 

min)  g 

BMAP-Montauk 
(long  IR) 
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X 

10  • 

-0.94 
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X 
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-  1.46 
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2.46 

X 
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-2.75 

X 
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X 
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X 
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X 
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X 
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X 
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X 
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X 

10 
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X 

10  ^ 
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X 

10  ' 
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3.01 

X 
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-2.34 

X 

10  ' 

5.35 

X 

10  ’  2 
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3.32 

X 

10 
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X 

10 
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X 
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ERIM  S23CH1 

1.21 

X 
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-8.08 

X 
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X 
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Figure  17  shows  each  of  the  images  in  Fig.  1  after 
filtering  was  applied.  To  see  the  results  of  the  process¬ 
ing  clearly,  only  a  portion  of  each  image  is  displayed. 
A  pi.xel  area  of  about  300  x  300  is  displayed  for  each 
image.  For  the  AFGL  scene  this  is  the  entire  image. 

The  images  were  created  by  first  taking  the  abso¬ 
lute  value  of  the  data.  Each  of  the  images  was  then 
scaled  from  0  to  255  with  0  equal  to  0.0  W/cm  Vsr, 
and  255  equal  to  the  maximum  filtered  radiance  in  the 
image.  Each  image  was  then  thresholded  to  be  con¬ 
sistent  with  succeeding  processing.  All  levels  below  the 
threshold  were  set  to  black;  all  levels  above  the  thresh¬ 
old  were  set  to  white.  The  thresholds  for  the  seven 
scenes  are  shown  in  Table  10  in  terms  of  the  noise 
equivalent  irradiance.  This  thresholding  will  be  dis¬ 
cussed  later. 

Using  this  thresholding,  any  negative  or  positive  out¬ 
puts  of  the  filter  that  are  greater  than  the  threshold 
are  white  in  the  images.  The  BMAP  image  in  Fig.  17a 
shows  that  about  50%  of  the  figure  is  greater  than  the 
threshold.  Figure  17b  is  mostly  black  with  a  much 
smaller  percentage  of  pixels  greater  than  the  thresh¬ 
old.  Figure  17c,  d,  and  e  shows  the  morning,  after¬ 
noon,  and  evening  scenes  of  Camp  A.P.  Hill,  after 
filtering.  Figure  1 7c  and  d  shows  only  a  small  percent¬ 
age  of  pixels  above  the  threshold.  Only  a  very  few  pix¬ 
els  are  above  the  threshold  in  Fig.  17e. 

One  bad  line  is  visible  in  the  Port  Hueneme  scene 
of  Fig.  17f.  The  high-pass  filtering  and  thresholding 
have  accented  this  bad  line  since  it  was  not  visible  in 
Fig.  If.  Also  shown  in  Fig.  17f  are  strong  edges  from 
the  buildings  in  the  scene.  The  ridges  in  the  mountains 
of  Fig.  17g  also  give  strong  edges  after  high-pass 
filtering. 

4.4  Filtering  and  Pulse  Discrimination  Performance 
Analysis 

After  applying  the  filters,  the  next  step  in  the  pro¬ 
cedure  is  a  pulse  discrimination  algorithm.  This  step 
adds  an  additional  false-alarm-rejection  capability  to 
suppress  sharp  but  extended  features  such  as  edges  and 
boundaries.' 

Figure  18  shows  the  pulse  width  discrimination  al¬ 
gorithm  for  the  BMAP  images.  The  lower  part  of  the 
figure  shows  a  point  targe  (a  near  impulse).  After 
filtering,  this  target  produces  two  pulses,  one  positive 
and  one  negative,  which  is  the  impulse  response  of  the 
bandpass  filter.  The  magnitude  of  the  spikes  and  the 
delay  between  the  maximum  and  minimum  are  func- 

■'I  .  M.  Howser.  "Investigation  of  Infrared  Background  flui 
ter."  in  fiesean  h  and  Development  Programs  at  The  John\  Hop¬ 
kins  University  Applied  Physics  Lahoratorv.  JHl  API 
yuarierly  Report  RQR'85-2  (Apr  Jun  I9KS). 


tions  of  the  order  and  the  cutoff  frequency  of  the  filter. 

The  upper  part  of  Fig.  18  shows  a  step  input.  This 
is  used  to  simulate  the  edge  of  a  cloud  as  seen  by  the 
sensor.  After  filtering,  this  artifact  also  produces  a 
positive  pulse  but  no  corresponding  negative  pulse. 
Similarly,  a  negative  step  would  produce  a  negative 
pulse  but  not  a  corresponding  positive  pulse. 

Any  post-filtering  threshold  set  to  detect  point  tar¬ 
gets  could  not  distinguish  between  the  point  target  and 
the  step  artifact  impulses.  To  perform  an  additional 
discrimination  between  point  and  extended  features, 
a  pulse-width  measurement  algorithm  is  used.  The 
pulse-width  discrimination  algorithm  searches  for  any 
local  maximum  above  a  specified  threshold.  (The  al¬ 
gorithm  actually  searches  for  a  maximum  or  minimum 
since  a  target  could  be  hotter  or  colder  than  the  sur¬ 
rounding  background.  We  will  focus  on  the  maximum 
in  this  discussion.)  After  finding  a  maximum,  succeed¬ 
ing  filter  output  is  searched  to  find  a  corresponding 
minimum  based  on  a  given  threshold.  If  a  maximum 
and  a  corresponding  minimum  are  both  found,  the 
point  is  declared  to  be  a  target;  therefore,  in  Fig.  18 
the  point  target  would  pass  the  discrimination  logic 
and  the  step  artifact  would  not  pass. 

Three  parameters  are  needed  for  the  pulse- 
discrimination  algorithm.  One  parameter  is  the 
minimum-to-maximum  ratio.  This  parameter  requires 
that  the  absolute  value  of  the  minimum  pulse  that  fol¬ 
lows  the  local  maximum  passing  the  detecting  thresh¬ 
old  must  be  of  some  minimum  size.  The  second 
parameter  is  the  delay  between  the  maximum  and  mini¬ 
mum  values.  Since  this  algorithm  searches  for  point 
targets,  the  delay  is  small.  The  final  parameter  that 
must  be  specified  is  the  detection  threshold.  To  avoid 
processing  sensor  noise  as  possible  targets,  any  pixel 
less  than  the  detection  threshold  is  not  processed. 

As  discussed  earlier,  the  selected  filter  determines 
the  theoretical  maximum  for  an  impulse  response  and 
the  delay  between  the  maximum  and  minimum.  For 
this  reason  the  thresholds  and  delays  used  in  this  anal¬ 
ysis  varied  depending  on  the  filters  used.  In  order  to 
determine  the  appropriate  thresholds  and  delays,  var¬ 
ious  inputs  were  processed  through  the  filters.  The  re¬ 
sponse  of  the  filters  to  these  inputs  was  used  to 
determine  the  thresholds  and  delays.  Since  each  set 
of  data  (BMAP,  AFGL,  and  ERIM)  had  different 
filters  this  selection  had  to  be  performed  for  each  data 
set. 

Since  the  BMAP  data  were  oversampled  by  the  seek¬ 
er,  simply  putting  an  impulse  into  the  filtering  would 
not  be  correct.  Before  filtering,  the  impulse  was  con¬ 
volved  with  the  point-source  characteristic  of  the  sen¬ 
sor,  spreading  the  input  into  more  than  one  output 
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Figur*  17  Images  after  filter  and  threshold,  (a)  BMAP  middle-wave  image  from  Bedford  after  bandpass 
filter  and  1 .5  NEI  threshold,  (b)  AFGL  scene  109  of  Monument  Valley  after  high-pass  filter  and  1 .5  NEI  thresh¬ 
old.  (c)  Morning  (scene  S01CH2).  (d)  Afternoon  (scene  S02CH2).  (e)  Evening  (scene  S03CH2).  (f)  ERIM  im¬ 
age  of  Port  Hueneme  (4. 5-5.5  ^m  spectral  band)  after  high-pass  filter  and  10.0  NEI  threshold  (scene  S23CH3). 
(g)  ERIM  image  of  mountains  by  Nellis  AFB  (3.9-4.7  >im  spectral  band)  after  high-pass  filter  and  3.0  NEI 
threshold  (scene  F8CH1). 
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Table  10  Thresholds  used  in  Fig.  17. 


Maximum  Threshold  Threshold 
Image  (W/cmVsr)  (W/cmVsr)  (LUT)  x  NEl' 


BMAP,  Bedford,  mid  IR 

1.18 

X 

10  ’ 

1.12 

X 

10'^ 

2 

1.5 

AFGL  {I09R) 

2.67 

X 

10  •’ 

5.60 

X 

10 

53 

1.5 

ERIMI  (S01CH2) 

3.96 

X 

10  - 

3.36 

X 

10  ^ 

22 

3 

(S02CH2) 

4.22 

X 

10  ' 

3.36 

X 

10“* 

20 

3 

(S03CH2) 

3.01 

X 

10  * 

3.36 

X 

10“^ 

29 

3 

ER1M3  (S23CH3) 

1.55 

X 

10  •* 

1.25 

X 

10  - 

21 

10 

ER1M5  (F8CH1) 

2.88 

X 

10  - 

1.63 

X 

10  ^ 

14 

3 

X  NEl  =  times  noise  equivalent  irradiance 


Scene 

Ma«  *  t.OO- 


Artifact  » . 


Time 


BMAP 

Bandpass  filter 


Threshold  logic 


Max  *  0>20S 


Fail 


Time 


Pass 


Figure  18  Pulse-width-discrimination  algorithm  for  BMAP  data. 


point  (performing  a  low-pass  function  on  the  impulse). 
The  result  of  the  convolution  was  filtered  with  the 
bandpass  filter  selected  for  the  BMAP  data.  In  addi¬ 
tion  to  the  impulse,  other  inputs  were  processed  in  the 
same  way.  After  the  inputs  were  processed,  informa¬ 


tion  on  the  threshold  and  delays  was  obtained.  Table 
1 1  shows  the  results  of  this  process. 

Table  1 1  shows  than  an  impulse  that  has  been  con¬ 
volved  with  the  BMAP  point-source  characteristic  has 
a  magnitude  of  about  0.25  times  its  original  value  af- 
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Table  11  Results  of  processing  impulses  with  BMAP  filters. 


Input 

signal 

JL~ 

JR_ 


I  (  '  J 

!  I  I  1 


I  I  I  I  I 
1  I  I  I  I 


JTh. 


Input  max 
(convolved) 

Output 

max 

Output 

min 

Max  -  min 

Delay 

Output  min 
Output  max 

255.0 

64.3 

-  51.0 

115.3 

3 

0.793 

468.7 

125.9 

-  91.2 

217.0 

4 

0.724 

681.4 

158.9 

-  125.5 

284.4 

4 

0.790 

793.6 

186.7 

-  145.7 

332.4 

4 

0.780 

881.3 

197.0 

-151.3 

348.4 

5 

0.768 

921.6 

200.0 

-  150.7 

350.7 

6 

0.754 

361.8 

93.8 

-  71.2 

165.9 

4 

0.759 

362.3 

95.2 

-  68.2 

163.4 

3 

0.716 

469.0 

110.4 

-  88.4 

198.8 

4 

0.801 

575.4 

139.8 

-  105.5 

245.3 

3 

0.755 

631.3 

156.4 

-118.6 

274.9 

4 

0.758 

574.9 

142.4 

-  108.4 

250.9 

4 

0.761 

ter  filtering  with  the  Butterworth  bandpass  filter.  Also 
the  ratio  of  the  output  minimum  to  maximum  is  0.793. 
There  is  a  3-sample  delay  between  the  output  maxi¬ 
mum  and  minimum  for  the  impulse  input. 

The  second  line  of  Table  1 1  shows  the  results  for 
a  two-pixel-wide  input.  For  this  input,  the  ratio  of  out¬ 
put  minimum  to  maximum  is  0.724  and  the  delay  be¬ 
tween  the  two  is  4  samples. 

The  logic  for  the  pulse  discrimination  algorithm  was 
determined  using  the  results  in  Table  1 1.  A  value  of 
0.70  was  chosen  for  the  ratio  of  the  pulse-measuring 
threshold  to  the  initial  detected  pulse.  This  value  means 
that  if  a  local  maximum  is  found  and  the  correspond¬ 
ing  minimum  is  less  than  0.70  times  the  maximum, 
that  point  will  not  be  chosen  as  a  target. 

From  Table  1 1  a  delay  of  3  or  4  samples  was  cho¬ 
sen  for  the  algorithm.  Any  local  maximum  and  mini¬ 
mum  that  is  found  that  has  a  delay  less  than  3  or 
greater  than  4  will  not  be  chosen  as  a  target.  This  re¬ 
striction  will  eliminate  wide  artifacts  (like  the  5-sample¬ 
wide  pulse  in  the  table). 


With  the  restrictions  described,  all  of  the  wide  fea¬ 
tures  in  Table  1 1  will  not  be  chosen  as  targets.  Only 
inputs  less  than  4  samples  wide  will  pass  this  logic. 

The  only  parameter  that  still  needs  to  be  defined 
for  the  BMAP  images  is  the  noise  threshold.  The 
threshold  of  images  was  set  depending  on  the  NEI  of 
the  seeker.  The  NEI  of  the  system  is,  by  definition, 
the  irradiance  required  for  a  signal  to  have  a  magni¬ 
tude  equal  to  the  rms  noise  within  the  system.  The  NEI 
of  the  BMAP,  AFGL,  and  ERIM  sensors  are  given 
in  Table  3. 

Since  the  threshold  is  defined  in  terms  of  the  stan¬ 
dard  deviation  of  the  noise  of  the  system,  the  percent¬ 
age  of  the  image  expected  to  be  greater  than  the 
threshold  can  be  computed  using  statistics  tables.  (This 
calculation  assumes  a  normal  distribution  of  the  out¬ 
put  levels.)  If  the  threshold  is  set  to  be  3  NEI,  0.26% 
of  the  pixels  are  expected  to  be  greater  than  the  thresh¬ 
old.  With  a  threshold  of  1.5  NEI,  13.4%  of  the  pix¬ 
els  would  be  expected  to  be  greater  than  the  threshold. 
For  the  analysis  below,  the  detection  threshold  was 
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typically  set  at  a  low  multiple  of  NEI  to  develop  de¬ 
tection  statistics  in  the  small  scenes  that  are  available. 

Figure  19  demonstrates  the  effect  of  pulse  discrimi¬ 
nation  on  one  line  of  the  Montauk  long-IR  image.  Fig¬ 
ure  19a  shows  the  first  line  of  the  first  frame  of  the 
long-IR  Montauk  image  before  any  filtering.  Figure 
19b  shows  the  same  line  of  the  image  after  the  band¬ 
pass  filter  was  applied.  Notice  that  both  axes  have 


0  150  300  450 

Pixel 


Figure  19  Line  1  of  BMAP  Montauk  Point  long- 
wavelength  background  data,  (a)  Unprocessed  line  of 
scene,  (b)  Line  of  scene  after  high-pass  filter  (with  ±  1 .5 
NEI  threshold),  (c)  Line  of  scene  after  high-pass  filter 
and  pulse  discrimination  (with  ±1.5  NEI  threshold). 


changed.  They-axis  changed  because  the  mean  of  the 
image  was  removed  by  the  filtering.  The  x-axis  changed 
because  the  image  had  to  be  padded  with  zeros  to  per¬ 
form  some  of  the  processing.  The  horizontal  lines  in 
Fig,  19b  show  the  ±  1 .5  NEI  threshold  applied  to  the 
data,  demonstrating  that  most  of  the  pixels  in  the  line 
are  excluded  from  processing  with  a  simple  threshold. 
The  pulse  discrimination  process  is  then  needed  to 
eliminate  the  remaining  pixels.  Figure  19c  shows  the 
same  line  of  the  image  after  pulse  discrimination;  all 
but  one  of  the  pulses  that  crossed  the  threshold  have 
been  eliminated  by  using  the  pulse  discrimination. 
Also,  the  remaining  pulse  barely  exceeds  the  thresh¬ 
old.  If  the  threshold  were  increased  to  2  x  NEI  no  pix¬ 
els  in  this  line  would  pass  the  discrimination  logic. 

No  convolution  of  an  input  pulse  with  the  sensor 
response  function  was  needed  for  the  AFGL  data  be¬ 
cause  no  oversampling  occurred  with  this  seeker.  Ex¬ 
cluding  this  step,  the  same  procedure  was  applied  to 
the  AFGL  data  as  was  applied  to  the  BMAP  data.  Var¬ 
ious  inputs  were  processed  with  the  AFGL  high-pass 
filter  and  their  results  were  tabulated.  Table  12  shows 
the  results. 

Table  12  shows  that  the  output  minimum-to- 
maximum  ratio  was  0.586  for  the  impulse  input.  The 
delay  between  the  maximum  and  minimum  was  1  pixel 
(  =  one  sample).  These  results  show  there  is  a  large 
difference  between  the  more  realistic  bandpass  filter 
impulse  of  the  oversampled  BMAP  data  and  the  high- 
pass  filtered  impulse  of  data  sampled  at  the  Nyquist 
rale  (i.e.,  AFGL  and  ERIM). 

The  parameters  for  the  pulse  discrimination  al¬ 
gorithm  were  calculated  from  the  results  in  Table  12. 
A  ratio  of  0.5  was  selected  for  the  minimum-to- 
maximum  pulse  height.  A  delay  of  1  or  2  pixels  was 
allowed  between  the  maximum  and  minimum.  This 
selection  of  parameters  only  passes  I  -  and  2-pixel-wide 
targets.  No  other  features  will  pass  this  logic. 

The  threshold  for  the  AFGL  data  was  not  as  easy 
to  pick  as  it  was  for  the  BMAP  data.  The  scaling  of 
the  images  was  not  set  to  the  NEI  of  the  system.  The 
PSD  of  each  image  was  needed  to  compute  the  NEI 
of  the  system .  The  PSD  of  each  image  flattened  at  the 
noise  level  of  the  system.  The  minimum  value  of  each 
PSD  in  the  flattened  region  was  obtained  from  the 
plots.  The  threshold  for  the  AFGL  images  was  set  at 
3.0  NEI.  Knowing  the  scaling  of  each  of  the  images, 
the  threshold  in  units  of  output  levels  was  calculated. 
Table  13  lists  the  thresholds  used  for  the  AFGL 
images. 

The  filtering  process  was  again  repeated  for  the 
ERIM  data.  Two  filters  had  to  be  applied  to  these  im¬ 
ages  because  the  ERIM  data  had  two  spatial  resolu- 
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Table  12  Results  of  processing  impulses  with  AFGL  filter. 


Input 

signal 

Input 

max 

Output 

max 

Output 

min 

Max  -  min  Delay 

Output  min 
Output  max 

J 

n_ 

256.0 

181.0 

-  106.0 

287.0 

1 

0.586 

STL 

256.0 

181.0 

-  150.0 

331.0 

2 

0.829 

J 

rm 

256.0 

181.0 

-  168.2 

349.2 

3 

0.928 

J 

RTTL 

256.0 

181.0 

-  175.7 

356.7 

4 

0.971 

J 

1  1  1  I  ' 

<  1  1  1 

1  1  1  1 

L  256.0 

181.0 

-  178.8 

359.8 

5 

0.988 

Table  13 

Thresholds  of  AFGL  images  for  pulse  discrimination. 

Scene 

NEI 

(W'/cm’) 

ME 

(W/em’/sr) 

Threshold 
(3  X  NEI) 
(output  levels) 

101 R 

5.8 

X  10 

8.4  X 

10  ^ 

13.35 

104R 

1.4 

X  10  " 

2.1  X 

10  • 

12.84 

105R 

4.9 

X  10  ’■ 

8.4  X 

10 

11.28 

I09R 

5.8 

X  10  " 

1.1  X 

10  ■ 

10.59 

lions.  Again,  \arious  inputs  were  processed  through 
the  filters  and  the  results  tabulated.  Table  14  shows 
the  results  of  th'  filtering.  The  first  section  of  this  ta¬ 
ble  is  for  the  ERIMl  data  (filtered  with  a  cut-on  of 
0. 10);  the  second  section  is  for  the  ERIM3  and  ERIM5 
data  (filtered  with  a  cut-on  of  0.070). 

The  results  in  Table  14  were  used  to  determine  the 
parameters  of  the  pulse-di.scrin. ’nation  logic  on  the 
ERIM  data.  The  ratios  for  ERIM3  images  were  slightly 
lower  than  the  corresponding  ratios  for  the  ERIMl 
images.  For  the  impulse  input,  the  ratio  of  the  out¬ 
put  minimum  to  output  maximum  was  0.84  for  the 
ERIMl  filters  and  0.75  for  the  ER1M3  and  ER1M5 
filters.  Looking  at  Table  14  it  can  be  seen  that  limit¬ 
ing  the  delay  between  maximum  and  minimum  to  I 
or  2  pixels  eliminates  some  of  the  wider  features  from 
being  considered  a  target.  It  would  be  desirable  to  set 
the  minimum-to-maximum  ratio  to  0.70  (just  below 


the  impulse  ratio);  however,  this  high  threshold  would 
eliminate  the  2-pixel-wide  pulse  in  line  5  of  both  sec¬ 
tions  of  the  table.  The  minimum-to-maximum  ratio 
was  therefore  set  to  0.50  in  order  to  avoid  eliminat¬ 
ing  possible  targets  from  consideration. 

The  last  parameter  to  be  specified  is  the  threshold 
for  both  sets  of  the  ERIM  data.  Two  methods  were 
used  in  computing  the  NEI  for  these  images.  The  first 
method  used  was  the  same  as  the  AFGL  calculation 
of  the  NEI.  The  minimum  of  the  PSD  was  assumed 
to  be  the  noise  power.  The  NEI  was  computed  from 
the  power.  The  second  method  used  the  noise  equiva¬ 
lent  differential  temperature  of  the  sensor.  Some  of 
the  ERIM  spectral  bands  had  a  noise  equivalent 
differential  temperature  (NEAT)  already  computed. 
Figure  20  shows  a  plot  of  the  NEAT  of  some  of  the 
bands  of  ERIM  data  versus  temperature.  An  average 
temperature  of  285  K  was  assumed  for  the  ERIM  im- 


Noise  equivalent  differential  temperature 
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Table  14  Results  of  processing  impulses  with  ERIM  filter. 
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Table  IS  ERIM  NEIs  and  thresholds. 


Scene 

NEI 

(W/cm^) 

FOV 

(circular) 

(mrad) 

3  X  NEI/(FOV) 
(W/cmVsr) 

SOI  CHI 

2.0 

3.68 

X 

10* 

(from  PSD) 

S01CH2 

4.47  X 

10’'^ 

2.0 

3.36 

X 

10* 

S01CH3 

2.0 

1.84 

X 

10’ 

S02CH1 

2.0 

3.68 

X 

10'* 

(from  PSD) 

S02CH2 

4.47  X 

10"'^ 

2.0 

3.36 

X 

10* 

S02CH3 

2.0 

1.16 

X 

10’ 

(from  PSD) 

S03CH2 

4.47  X 

10'^ 

2.0 

3.36 

X 

10'* 

S03CH3 

2.0 

1.50 

X 

10“’ 

(from  PSD) 

S23CH1 

2.5 

5.36 

X 

10* 

(from  PSD) 

S23CH3 

7.80  X 

10 

2.5 

3.74 

X 

10'* 

S23CH4 

4.63  X 

10"“ 

2.5 

2.20 

X 

10*’ 

F8CHI 

3.40  X 

lO"'^ 

2.5 

1.63 

X 

10'* 

F8CH2 

8.83  X 

10 

2.5 

4.24 

X 

10* 

F8CH3 

4.58  X 

10" 

2.5 

2.20 

X 

10’ 

Table  16  Pulse  discrimination  parameters. 


Data 

source 

Output  min 
Output  max 

Delay 

(pixels) 

Threshold 
(X  NEI) 

BMAP 

0.70 

3  or  4 

1.5 

AFGL 

0.50 

1  or  2 

3.0 

ERIM 

0.50 

1  or  2 

3.0' 

'Threshold  started  at  3.0  but  had  to  be  increased 
for  some  images. 


in  Table  16.  The  results  of  the  processing  are  listed 
in  Table  17,  which  lists  the  maximum  and  minimum 
of  each  of  the  images  after  this  processing.  Also  list¬ 


ed  in  Table  17  is  the  difference  between  the  maximum 
and  minimum.  The  threshold  used  (in  W/cm^/sr)  is 
listed  in  the  following  column.  The  threshold  in  NEI 
is  listed  in  parentheses  in  the  same  column.  The  per¬ 
centage  of  samples  that  passed  this  algorithm  is  listed 
in  the  last  column. 

Some  of  the  images  were  processed  twice  through 
the  pulse-discrimination  algorithm.  If  too  many  pix¬ 
els  passed  the  logic  the  first  time,  the  threshold  was 
raised  and  the  image  was  reprocessed.  The  results  for 
all  of  these  runs  are  also  listed  in  Table  17. 

Figure  21  shows  the  seven  sample  images  after  filter¬ 
ing  and  pulse  discrimination.  The  same  procedure  as 
was  used  for  the  images  in  Fig.  17  was  followed  to 
threshold  these  images.  Also,  the  same  portion  of  each 
image  was  displayed  to  make  possible  direct  compar¬ 
isons  between  the  images. 
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Table  17  Statistics  of  data  after  filtering  and  pulse  discrimination  in  units  of  W/cm^/sr. 


Image 

Max 

Min 

(max 

A 

- 

min) 

Threshold 
(X  NEI) 

Pixels 

passed 

C^o) 

BMAP-Montauk 

3.68 

X 

10’* 

-6.38 

X 

10* 

1.01 

X 

10 

5 

1.38 

X 

10 

"  (1.5) 

0.10 

(Far  IR) 

BMAP-Montauk 

1.23 

X 

10"‘ 

-6.0 

X 

10  ’ 

1.83 

X 

10 

6 

2.75 

X 

10 

’  (1.5) 

0.10 

(Mid  IR) 

BMAP-Bedford 

6.18 

X 

10  '' 

-5.44 

X 

10 

1.16 

X 

10 

5 

1.5 

X 

10 

(1.5) 

0.04 

(Far  IR) 

BMAP-Bedford 

1.25 

X 

10  * 

-  1.27 

X 

10  * 

2.52 

X 

10 

6 

8.38 

X 

10 

*  (1.5) 

0.18 

(Mid  IR) 

AFGL  101 R 

1.34 

X 

10 

-  1.19 

X 

10  " 

2.52 

X 

10 

4 

5.61 

X 

10 

’  (1.5) 

1.02 

1.34 

X 

10  ^ 

-  1.19 

X 

10  " 

2.52 

X 

10 

4 

1.12 

X 

10 

"  (3.0) 

0.01 

AFGL  I04R 

1.80 

X 

10  ■' 

-2.43 

X 

10  " 

6.78 

X 

10 

4 

1.35 

X 

10 

^  (1.5) 

0.07 

0 

0 

0 

2.71 

X 

10 

^  (3.0) 

0 

AFGL  I05R 

1.69 

X 

10  ^ 

-2.77 

X 

10  " 

4.46 

X 

10 

4 

4.47 

X 

10 

’  (1.5) 

2.02 

1.69 

X 

10" 

-2.77 

X 

10  " 

4.46 

X 

10 

4 

0.48 

X 

10 

’  (3.0) 

0.05 

AFGL  I09R 

1.37 

X 

10  " 

-  1.29 

X 

10  " 
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The  image  in  Fig.  21a  shows  that  a  very  small  per¬ 
centage  of  the  points  passed  the  pulse-discrimination 
logic.  The  AFGL  image  (Fig.  21b)  does  not  show  a 
large  reduction  in  the  number  of  detected  pixels;  how¬ 
ever,  a  reduction  in  the  number  of  detected  pixels  is 
still  evident.  Figure  21c  and  d  shows  a  large  decrease 
in  the  number  of  pixels  above  the  threshold.  Only  a 
very  few  pixels  in  Fig.  21e  passed  the  threshold. 

Figure  21f  shows  the  same  area  as  Fig.  17f.  Only 
one  or  two  of  the  building  outlines  passed  the  pulse 
discrimination.  Many  of  the  smaller  points  also  passed 
the  logic.  The  bad  scan  line  is  still  visible  in  this  figure. 

Finally,  Fig.  21g  shows  the  mountains  in  the  ERIM5 
image.  A  comparison  with  Fig.  17g  shows  that  some 
of  the  mountain  ridges  have  been  removed  from  the 
image.  However,  many  of  the  single  pulses  were  not 
removed  from  the  image,  a  condition  that  would  cause 
many  false  alarms  with  this  algorithm. 

4.5  Cross-Scan  Comparison 

The  pulse-discrimination  algorithm  discussed  in  the 
previous  section  reduces  the  required  target  irradiance 
by  analyzing  the  pulse  shapes  along  the  scan  direction. 
A  cross-scan  comparison  checks  neighboring  scan  lines 
to  find  features  that  extend  across  several  scan  lines. 
If  this  analysis  were  performed,  the  number  of  false 
alarms  could  be  reduced  further. 

A  cross-scan  comparison  is  more  difficult  to  per¬ 
form  than  the  pulse-discrimination  algorithm  because 
the  extended  features  could  be  at  any  angle  relative 
to  the  scan  direction.  Only  checking  pixels  perpendic¬ 
ular  to  the  scan  direction  would  not  be  sufficient.  A 
method  to  recognize  extended  features  at  any  angle 
would  need  to  be  developed. 

No  cross-scan  algorithm  was  implemented  for  this 
analysis.  The  thresholded  images  in  Fig.  21  were 
checked  visually  for  any  extended  features.  The  in¬ 
spection  can  tell  qualitatively  whether  the  required  tar¬ 
get  irradiance  would  be  reduced  further  by  applying 
a  cross-scan  comparison  after  the  pulse-discrimination 
algorithm. 

For  Fig.  21a  through  e,  it  is  doubtful  that  a  cross¬ 
scan  comparison  would  further  reduce  the  number  of 
pixels  that  passed  the  discrimination  logic.  Few,  if  any, 
of  the  pixels  in  these  images  form  a  continuous  line. 
However,  in  the  port  and  mountain  scenes  in  Fig.  21  f 
and  g  the  cross-scan  comparison  could  eliminate  some 
of  the  false  alarms. 

In  Fig.  2 if,  one  edge  of  each  of  two  buildings  is 
roughly  outlined  by  false  alarms.  A  cross-scan  com¬ 
parison  could  find  these  outlines  and  eliminate  those 
samples  from  being  considered  possible  targets. 

A  cross-scan  comparison  would  also  be  useful  in 


Fig.  21g.  Many  short  vertical  lines  are  visible  in  this 
image.  Also,  a  long,  continuous  feature  is  found  on 
the  right  edge  of  the  image.  If  these  features  could  be 
removed  from  being  considered  possible  targets,  many 
of  the  false  alarms  in  this  image  would  be  removed. 

4.6  Comparison  of  Processing  Techniques 

The  filtering  and  pulse  discrimination  logic  were  ap¬ 
plied  to  these  images  to  reduce  the  noise  level  of  the 
background,  thereby  improving  the  signal-to-clutter 
ratio.  One  measure  of  the  improvement  is  the  mini¬ 
mum  detectable  signal.  For  the  unprocessed  images, 
the  smallest  signal  that  could  be  detected  above  the 
clutter  would  have  to  be  greater  than  the  maximum 
pixel  in  the  image  (or  less  than  the  minimum  pixel  if 
both  a  maximum  and  minimum  threshold  were  set). 

The  minimum  detectable  signal  for  the  images  pro¬ 
cessed  through  only  the  filtering  or  processed  through 
the  filtering  and  pulse  discrimination  can  be  calculat¬ 
ed  from  the  maximum  and  minimum  values  of  these 
images.  It  is  shown  in  Tables  11,  12,  and  14  that  the 
maximum  value  of  an  impulse  is  reduced  by  the  filter¬ 
ing.  For  example,  in  Table  1 1  the  maximum  of  the 
input  is  255.0  while  the  maximum  of  the  output  is  64.3. 
This  gives  an  output-to-input  ratio  of  about  0.25.  Be¬ 
cause  of  this  effect,  the  minimum  detectable  signal  for 
images  that  have  been  filtered  is  calculated  by  divid¬ 
ing  the  maximum  (and  minimum)  values  of  the  im¬ 
ages  by  the  ratio  of  maximum  output  to  maximum 
input.  This  calculation  produces  the  required  signal 
into  the  filtering  to  produce  the  given  maximum  or 
minimum  after  filtering.  Table  18  shows  the  factors 
used  to  calculate  the  minimum  detectable  signal.  These 
values  were  computed  from  the  impulse  results  in  Ta¬ 
bles  II,  12,  and  14. 

Table  19  lists  the  minimum  detectable  signal  for  each 
of  the  images  in  the  study  before  any  processing,  af- 


Table  18  Filtering  output-to-input 
ratios. 


Data 

source 

Output/input 

BMAP 

0.252 

AFGL 

0.707 

ERIMl 

0.578 

ER1M3,  ERIM5 

0.625 

THE  X>HNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


Table  19  Minimum  detectable  signal  in  units  of  W/cm^/sr. 


Before  processing 

After  filter 

After  filter 
and  pulse 
discrimination 

Image 

Max  -/i  Min  -fi 

Max  Min 

Max  Min 
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(Far  IR) 

1.1 

X 

10  '' 

-9.0 

X 

10-’ 
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10’ 
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X 

10’ 
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1  (Mid  IR) 

1.75 
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10’ 

-9.4 

X 

10* 
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X 

10"‘ 
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10"* 
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1.0 

X 

10-“ 
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10-“ 
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2.45 

X 

10"’ 

-2.16 

X 

10’ 

5  BMAP-Bedford 

1  (Mid  IR) 

8.7 

X 

10“’ 

-3.24 

X 

10"’ 

4.68 

X 

10"’ 

-1.71 

X 

10"’ 
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X 

10-* 

-5.04 

X 

10* 

{  AFGL  101 R 

1.0 

X 

10-’ 

-4.5 

X 

10-“ 

3.03 

X 

10“ 
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X 

10“ 

1.90 

X 

10-“ 

-1.68 

X 

10“ 
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6.2 

X 

10-“ 

-1.07 

X 

10"’ 

8.67 
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10"“ 

-3.44 
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10“ 
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X 

10“ 

-3.44 

X 

10“ 
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4.1 

X 
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X 
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-4.8 
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3.51 

X 

10“ 
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X 
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1.93 

X 

10-“ 

-1.82 

X 

10-“ 
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1 

1.99 

X 
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-2.67 

X 
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X 

10"’ 
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X 

10-’ 

4.62 

X 

10-’ 

-4.61 

X 

10"’ 

ERIM  F8CH2 

1.59 

X 
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-1.27 

X 

10"’ 

2.32 

X 
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-2.21 

X 
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1.65 

X 
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-2.21 

X 

10"’ 
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4.59 

X 
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-4.3 

X 
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2.98 

X 
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X 
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X 
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-2.40 

X 
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X 
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X 
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X 
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b 

X 
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X 
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X 
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f 
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X 
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1.8 

X 

10  ’ 
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-5.25 

X 

10  ’ 
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10  ’ 

-4.0 

X 

10  “ 

6.52 

X 

10  “ 

-6.68 
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X 
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9.0 

X 
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-4.0 

X 
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X 
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X 
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5.17 

X 

10  “ 
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X 

10  “ 
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2.22 

X 

10  “ 

-3.8 

X 

10  ’ 

1.94 

X 
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X 

10  “ 

1.74 

X 

10  “ 

-1.24 

X 

10  “ 

1  ERIM  S23CH3 

2.4 

X 

10  “ 

-6.8 

X 

10"’ 

2.48 

X 

10  “ 

-  1.5 

X 

10  “ 

1.68 

X 

10  “ 

-  1.46 

X 

10  “ 

\  ERIM  S23CH4 

2.2 

X 

10  ’ 

-7.5 

X 

10“ 

1.78 

X 

10  ’ 

-1.92 

X 

10  ’ 

1.78 

X 

10  ' 

-  1.92 

X 

10  ’ 
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ter  filtering,  and  after  filtering  and  pulse  discrimina¬ 
tion.  Since  the  means  of  the  images  are  removed  by 
the  filtering  the  means  of  the  before-processing  im¬ 
ages  have  been  subtracted  from  the  maximum  and 
minimum  values.  Therefore,  this  table  shows  the 
spread  about  the  mean  radiance  in  the  before-proces¬ 
sing  column  and  the  spread  about  0.0  in  the  remain¬ 
ing  columns. 

Figure  16  shows  the  relationship  between  the  before¬ 
processing  results  and  the  after-filtering  results.  Since 
the  low-frequency  features  are  present  in  the  original 
images,  any  spike  present  in  the  image  might  be  hid¬ 
den  by  the  low-frequency  features.  Therefore,  the  max¬ 
imum  and  minimum  values  in  this  column  might  not 
relate  to  spikes  in  the  images.  On  the  other  hand,  the 
low-frequency  features  are  removed  from  the  filtered 
images.  The  maximum  and  minimum  values  are  neces¬ 
sarily  from  spikes  in  the  image.  For  these  reasons  the 
maximum  or  minimum  after  filtering  could  be  great¬ 
er  than  in  the  original  image. 

Figure  22  displays  the  results  of  Table  19.  The  fig¬ 
ure  shows  only  the  positive  minimum  required  irradi- 


ance  values  for  the  BMAP,  AFGL,  and  ERIM  images. 
The  required  irradiance  before  processing,  after  filter¬ 
ing,  and  after  filtering  and  pulse  discrimination  are 
plotted  with  different  symbols.  Also  plotted  on  this 
figure  are  the  minimum  required  irradiance  values 
from  the  application  of  the  matched  filter  to  the  PSDs 
of  the  BMAP  and  AFGL  images. 

To  analyze  the  results  in  Table  19  and  Fig.  22  the 
data  must  first  be  broken  into  subsets.  The  BMAP  im¬ 
ages  were  all  of  cloud  backgrounds.  Table  19  shows 
that  the  minimum  detectable  signal  of  these  images 
was  reduced  by  the  filtering.  For  all  four  of  these  im¬ 
ages  the  pulse  discrimination  algorithm  further  reduced 
the  required  signal.  Figure  23  shows  a  histogram  of 
these  results.  Three  histograms  of  the  Bedford  mid- 
IR  data  are  presented  in  this  figure — the  original  im¬ 
age  (with  the  mean  subtracted),  the  image  after  filter¬ 
ing  with  a  first-order  high-pass  and  a  second-order 
low-pass  filter,  and  the  image  after  filtering  and  pulse 
discrimination.  The  figure  shows  a  large  reduction  in 
the  maximum  and  minimum  owing  to  the  filtering. 
The  smaller  reduction  owing  to  the  pulse  discrimina¬ 
tion  is  also  seen. 
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Radiance  (W/cm^/sr) 

Figure  23  Histograms  of  BMAP  mid-IR  Bedford 
scene. 


Figure  24a  and  b  further  demonstrates  the  effect  of 
the  pulse  discrimination  and  filtering  on  the  IR  scenes. 
This  figure  shows  a  cumulative  histogram  for  the  un¬ 
processed  image,  the  image  after  filtering,  and  after 
filtering  and  pulse  discrimination.  Figure  24a  shows 
these  histograms  for  AFGL  scene  I09R  and  Fig.  24b 
shows  the  histograms  for  ERIM  (S01CH2).  The  or¬ 
dinate  of  these  plots  is  the  cumulative  frequency  of 
occurrence.  Two  scales  are  drawn  along  the  abscissa — 
the  irradiance  and  the  NEI.  This  plot  indicates  where 
the  NEI  or  irradiance  threshold  must  be  set  to  get  no 
false  alarms  for  these  two  scenes.  These  plots  again 
show  the  decrease  in  the  required  irradiance  owing  to 
filtering  and  filtering  with  pulse  discrimination. 

Figure  22  shows  that  for  the  BMAP  scenes  the 
matched  filter  usually  requires  a  lower  target  irradi¬ 
ance  than  is  required  by  the  Butterworth  filter.  How¬ 
ever,  the  Butterworth  filter  followed  by  the  pulse- 
discrimination  algorithm  usually  performs  better  (i.e., 
requires  a  lower  target  irradiance)  than  the  matched 
filter.  This  unexpected  result  is  probably  attributable 
to  the  higher  threshold  used  in  the  matched  filter  anal¬ 
ysis.  The  matched-filter  threshold  was  much  higher 
than  needed  to  eliminate  all  false  alarms  in  these  small 
scenes. 

Scene  101 R  was  the  only  processed  AFGL  image  that 
contained  only  clouds.  The  results  in  Table  19  show 
that  the  filtering  and  pulse  discrimination  did  reduce 
the  minimum  detectable  signal.  The  remaining  AFGL 
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- Unprocessed  (with  mean  removed) 
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- Filtered  and  pulse  discrimination 

Figure  24  Cumulative  histograms  showing  reduction 
of  maximum  irradiance  and  NEI  from  filtering  and  pulse 
discrimination,  (a)  Histograms  of  AFGL  Monument  Val¬ 
ley  and  cloud  scene  (109).  (b)  Histograms  of  ERIM  Camp 
A.P.  Hill  morning  scene  4. 5-5.5  pm  spectral  band) 
(scene  S01CH2). 


scenes  had  clouds  and  terrain  in  the  images.  Scene  I09R 
shows  a  decrease  in  the  required  signal  for  the  filter¬ 
ing  and  the  pulse  discrimination.  Scenes  i04R  and  I05R 


46 


B 


•5 


a» 


S 


Si 


THE  X>HNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


show  a  reduction  owing  to  the  filtering,  but  the  pulse 
discrimination  does  not  further  reduce  the  negative  re¬ 
quired  signal.  This  result  means  that  there  were  large 
negative  pulses  in  the  images  that  passed  the  filtering 
and  the  pulse  discrimination.  This  result  is  expected 
to  be  more  common  for  terrain  images  than  for  cloud 
images. 

For  all  four  of  the  AFGL  scenes,  the  Butterworth 
filter  followed  by  the  pulse  discrimination  performed 
better  than  the  matched  filter  applied  to  the  PSD. 

The  ERIM  data  must  be  broken  into  many  parts 
for  analysis.  The  ERIMS  images  used  in  this  analysis 
IF8CH1,  F8CF12,  and  F8CH3)  were  of  mountains. 
The  results  in  Table  19  show  that  in  some  cases  the 
pulse  discrimination  further  reduced  the  minimum  re¬ 
quired  signal,  but  the  value  often  remained  almost  the 
same. 

The  set  of  images  from  ERIMl  (SOI  through  S03) 
can  be  discussed  by  band.  One  of  the  channel  1  im¬ 
ages  shows  very  little  decrease  in  the  minimum  detect¬ 
able  signal  when  pulse  discrimination  is  added. 
Channel  2  images,  on  the  other  hand,  all  show  a  sig¬ 
nificant  decrease  in  false  alarms  when  the  pulse  dis¬ 
crimination  is  added.  Some  of  the  channel  3  images 
show  a  decrease  in  the  maximum  or  minimum  clutter 
pulses  but  other  images  show  little  or  no  change  be¬ 
cause  of  the  pulse  discrimination. 

Scene  23  from  ERIM3  had  two  bad  lines  in  the  im¬ 
age.  Figure  25  shows  three  lines  of  the  image.  Figure 
25a  shows  a  normal  line  in  channel  3  (line  1)  of  the 
image.  Figure  25b  and  c  shows  lines  18  and  641  of 
channel  3.  These  plots  show  that  there  is  something 
wrong  with  these  two  lines  of  data,  a  problem  that 
was  not  found  until  after  the  pulse  discrimination  al¬ 
gorithm  was  applied  and  hundreds  of  points  from  each 
of  these  lines  passed  the  discrimination  logic.  Chan¬ 
nels  3  and  4  of  scene  23  had  bad  data  in  both  of  these 
lines.  Channel  1  had  bad  data  in  only  line  641.  Be¬ 
cause  of  these  bad  data,  the  results  in  Table  19  and 
Fig.  22  are  not  expected  to  show  a  significant  decrease 
resulting  from  the  pulse  discrimination. 

The  results  from  Table  19  can  be  related  to  the 
results  from  Table  17.  In  Table  19  the  BMAP  data 
(clouds)  showed  a  large  decrease  in  the  minimum  de¬ 
tectable  signal  owing  to  filtering  and  pulse  discrimi¬ 
nation.  The  results  in  Table  17  agree  with  these  results 
in  that  less  than  0.2%  of  the  pixels  in  the  images  passed 
the  discrimination  logic  with  the  threshold  set  at  1.5 
NEI.  Many  more  of  the  AFGL  pixels  passed  the  log¬ 
ic.  This  result  can  again  be  associated  with  the  scene 
content. 

Fewer  points  passed  the  pulse  discrimination  logic 
for  channel  2  of  ERIM  scenes  SOI,  S02,  and  S03  than 
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Figure  25  Lines  of  images  from  ERIM  Port  Hueneme 
scene  in  4,5-5.5  j^m  spectral  band  (scene  S23CH3) 
showing  bad  scan  lines,  (a)  Line  1  of  scene  with  no  bad 
data,  (b)  Line  18  of  scene  with  entire  line  of  bad  data, 
(c)  Line  641  of  scene  with  partial  line  of  bad  data. 


for  channels  1  and  3,  consistent  with  the  results  in  Ta¬ 
ble  19. 

The  bad  scan  lines  in  scene  S23  affect  the  results 
in  Table  17.  To  get  about  0.3%  of  the  points  to  pass 
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the  logic,  the  threshold  had  to  be  raised  to  20  NEl  for 
one  of  the  channels,  agreeing  with  the  results  in  Ta¬ 
ble  19  where  little  to  no  improvement  was  gained  with 
the  pulse  discrimination.  This  result  shows  that  most 
of  the  bad  points  passed  the  filtering  and  the  pulse 
discrimination. 

4.7  Implementation  of  the  Scanning  Sensor 
Processing 

The  previous  sections  of  this  leport  described  the 
applications  of  clutter  suppression  algorithms  to  dig¬ 
ital  infrared  images.  This  section  briefly  summarizes 
implementation  of  the  algorithms  in  a  missile  guid¬ 
ance  system. 

Figure  26  shows  a  block  diagram  of  one  potential 
seeker-processing  configuration  that  implements  the 
scanning-sensor  algorithms  developed  in  this  report. 
The  detected  scene  (plus  target)  is  preamplified  and 
prefiltered.  The  prefilter  suppresses  high  frequency 
(out  of  the  seeker  bandpass)  and  may  also  be  used  to 
suppress  the  detector  DC  output.  The  resultant  video 
signal  from  each  channel  is  then  digitized  in  an  A/D 
converter.  The  converter  has  a  dynamic  range  of  14 
to  16  bits  (or  as  few  as  8  bits  for  AC-coupled  long- 
wavelength  sensors),  and  samples  at  two  to  three  times 
the  Nyquist  rate.  The  output  digital  signal  is  then  avail¬ 
able  for  clutter  suppression  algorithm  implementation. 

The  digital  signal  is  first  processed  with  a  bandpass 
filter  that  is  matched  to  the  clutter  spectrum  and  the 
spatial  filter  at  the  scanning  detector.  Filtered  output 
is  then  threshoided  at  approximately  five  to  ten  times 
the  sensor  NEl  to  eliminate  most  of  the  internal  noise. 
Data  above  threshold  are  a  detection.  Data  below 


threshold  are  perserved  for  use  in  the  pulse-width  dis¬ 
crimination  algorithm.  In  general,  both  positive  and 
negative  thresholds  are  implemented  to  detect  targets 
of  both  contrasts. 

After  thresholding,  each  detection  is  pulse-width 
tested  by  looking  for  the  corresponding  pulse  of  the 
opposite  sign,  appropriate  magnitude,  and  appropri¬ 
ate  delay.  The  exact  implementation  of  this  pulse-width 
discriminator  is  dependent  on  the  filtering  and  noise 
level.  Tolerance  on  the  magnitude  of  the  second  pulse 
for  point-target  confirmation  can  be  computed  readily 
from  the  measured  first  pulse  magnitude  and  the 
known  (or  real-time  measured)  system  noise.  At  this 
point  the  data  are  compressed  and  only  detections  are 
preserved. 

The  final  step  in  the  discrimination  algorithm  is 
comparison  of  the  detected  output  at  adjacent  detec¬ 
tor  channels.  Detections  (passing  threshold  and  pulse- 
width  discrimination)  that  span  typically  three  or  more 
adjacent  channels  within  two  or  three  time  samples 
of  one  another  are  declared  false  targets.  If  the  de¬ 
tectors  are  offset  and  if  they  overlap  across  the  scan 
direction,  the  detection  of  false  targets  is  made  easi¬ 
er.  With  this  configuration,  a  false  target  would  be 
detected  in  three  or  more  contiguous  channels.  This 
would  eliminate  the  necessity  for  checking  channels 
two  or  three  time  samples  apart. 

Outputs  of  all  the  detector  channels  may  also  be 
used  to  raise  the  detection  threshold  adaptively,  on 
the  basis  of  occurrence  of  many  false  alarms  within 
a  scan  frame  (or  segment  thereof).  When  the  number 
of  detections  per  frame  becomes  large,  a  momentary 
increase  in  detection  threshold  may  be  necessary  to 
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Preamp  Prefilter  Digitizer 
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Digital  Detection  Pulse-width 
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Figure  26  Implementation  of  the  scanning  sensor  processing. 
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discriminate  against  clutter  patches  during  search.  to  avoid  the  raising  of  the  threshold  while  tracking 
Once  in  track,  different  algorithms  would  be  used  to  a  valid  target, 
discriminate  against  the  appearance  of  decoys  so  as 


5.0  CONCLUSIONS 


The  performance  of  a  scanning  sensor  in  back¬ 
ground  clutter  has  been  investigated.  A  number  of 
background  scenes,  supplied  by  three  separate  mea¬ 
surement  programs,  have  been  characterized  in  terms 
of  artifacts,  imagery,  scene  radiance  statistics,  and 
power  spectral  density. 

Available  scene  data  were  processed  in  several  ways 
to  determine  the  minimum  target  signature  that  can 
be  detected  in  particular  scenes.  Two  techniques  were 
used  to  derive  performance.  The  first  used  the  scene 
power  spectral  density  to  estimate  required  target  ir- 
radiance  as  a  function  of  sensor  resolution  (IFOV)  as¬ 
suming  a  fixed  threshold  based  on  a  matched-filter 
(maximum  signal-to-noise)  analysis. 

Since  the  matched-filter  results  are  based  on  assump¬ 
tions  that  are  not  met  by  real  scenes,  direct  process¬ 
ing  of  the  scenes  is  also  done.  The  analysis  includes 
three  levels  of  scene  processing  in  much  the  same  way 
as  a  sensor  signal  processor  detection  system.  These 
three  levels  are: 

1 .  A  high-pass,  first-order  digital  filter  that  whit¬ 
ens  the  clutter  power  spectral  density  and  approx¬ 
imately  maximizes  the  target  signal-to-clutter 
ratio.  The  result  is  directly  comparable  to  the 
matched-filter  results  but  gives  a  better  estimate 
of  false-alarm  rate.  The  filter  provides  good  dis¬ 


crimination  against  large,  low-frequency  scene 
features. 

2.  A  high-pass  filter,  as  above,  but  in  addition  there 
is  a  pulse-width  discriminator.  The  discrimina¬ 
tor  rejects  the  high-frequency  output  of  filtered 
edges  while  passing  the  high-frequency  output 
of  a  target  or  smaller  feature.  This  technique  is 
shown  to  reduce  significantly  the  false  alarms 
produced  by  a  scene. 

3.  A  high-pass  filter,  pulse-width  discriminator,  and 
cross-scan,  multiple  channel  logic.  Multiple- 
channel  logic  enhances  the  point-discrimination 
capability  of  the  system  by  determining  the  cross¬ 
scan  spatial  extent  of  linear  features  and  by  re¬ 
jecting  those  that  are  significantly  greater  in 
cross-scan  extent  than  is  expected  of  point  tar¬ 
gets.  This  feature  was  particularly  useful  in  im¬ 
proving  performance  against  certain  terrain 
scenes  that  had  significant  linear  feature  content. 

The  results  of  the  direct-processing  methods  are  illus¬ 
trated  in  images,  graphs,  and  tables  that  compare  per¬ 
formance  of  these  techniques  as  well  as  the  effects  of 
spectral  bands,  scene  types,  and  data  spatial  resolu¬ 
tion.  The  analysis  is  concluded  with  a  discussion  of 
the  implementation  of  the  background  processing  in 
a  missile  guidance  system. 
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Appendix  A 


BMAP  Data  Description 


Additional  information  on  the  Navy  Backgrounds 
Measurement  and  Analysis  Program  (BMAP)  data  are 
contained  in  this  appendix.  Included  is  a  list  of  the 
characteristics  of  the  sensor,  the  long-  and  mid-wave 
spectral  responses,  and  curves  showing  the  long-  and 
mid-wave  point-source  characteristics. 

The  objective  of  the  BMAP  program  is  to  obtain 
clutter  data  for  design  of  air-  and  surface-based  IR 
search  and  track  sensors.  The  data  obtained  by  APL 
were  the  initial  surface-based  measurements  obtained 
with  a  Raytheon-owned  instrument.  Currently,  the 
described  instrument  is  collecting  data  from  a  P-3  air¬ 
craft.  An  advanced  instrument,  with  significantly  more 
detector  channels,  smaller  IFOV,  and  a  more  capa¬ 
ble  data  acquisition  system  is  being  developed  to  sup¬ 
port  future  measurements. 

Primary  characteristics  of  the  BMAP  sensor  of  in¬ 
terest  to  this  study  are  the  high-sensitivity,  simulta¬ 
neous  dual-band  capability,  and  the  availability  of 
oversampled  data.  Table  A-1  gives  the  basic  charac¬ 
teristics  of  the  BMAP  sensor.  Figure  A-l  gives  the 
spectral  response  of  the  mid-IR  and  long-IR  spectral 
bands.  Figure  A-2  gives  the  point-source  response  of 
typical  mid-IR  and  far-IR  detector  channels. 


Two  data  sets  were  available.  The  first  measurement 
series  was  conducted  at  Montauk  Point,  Long  Island, 
on  August  9-19,  1983.  These  data  are  primarily  over¬ 
cast  sky  and  have  less  contrast  than  the  later  Bedford 
data.  Fifty  frames  of  data  taken  in  the  8/8  (8  lines 
of  mid-IR  and  the  corresponding  8  lines  of  far-IR  data) 
format  on  August  10,  1983  were  provided.  The  sec¬ 
ond  set  of  data  was  taken  on  September  13,  1984  at 
Bedford,  Massachusetts.  Seven  scenes  of  a  hole  in  a 
cloud  deck  were  taken.  Table  A-2  summarizes  the 
BMAP  data. 

Calibration  constants  to  convert  digital  data  counts 
to  radiance  are  as  follows: 

Montauk  Point 

Long-IR  ~9.2  x  10  ’  W/cm^/sr  per  count 
Mid-IR  ~  1.8  X  10  ^^  W/cm^/sr  per  count 


Bedford 

Long-IR 

Mid-IR 


1.000  X  10  W/cmVsr  per  count 
5.587  X  10  "*  W/cm*/sr  per  count 


Multiplying  the  recorded  data  by  these  factors  recovers 
the  measured  radiance. 


v.v. 
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Table  A-1  Characteristics  of  the  BMAP  dual-band  IR  measurement  system. 


Parameter 

Value 

Collecting  area  (3. 9-4. 8  /xm)  (7.6-11.3  /xm) 

80  cm^ 

Optical  transmission  (3.9-4. 8  nm) 

(7.6-11.3  ;xm) 

85% 

80% 

Detector  collecting  efficiency  (3. 9-4. 8  /xm) 

(7.6-11.3  /xm) 

80% 

Total  optical  efficiency  (3. 9-4. 8  /xm) 

(7.6-11.3  /xm) 

68% 

64% 

Detector  area  (3. 9-4. 8  /xm) 

2.58  X  10"^  cm^ 

Detector  DJv  (1  Hz  -  1  kHz)  (3. 9-4. 8  /xm) 

(1  Hz  -  1  kHz)  (7.6-11.3  /xm) 

2  X  10"  cm  Hz'^VW 

1.5  X  10'“  cm  Hz'VW 

Electronic  3  dB  bandwidth  (3.9-4. 8  /xm) 

(7.6-11.3  /xm) 

0.0  Hz  -  1  kHz 

0.5  Hz  -  1  kHz 

Number  of  channels  single  color 

Two-color  (per  color) 

16"> 

8 

Total  FOV'  elevation 
azimuth 

5.33  mrad  or  0.31° 

2.8° 

Instantaneous  FOV 

0.33  X  0.33  mrad 

NEI  (3.9-4.8  /xm)  (PV  InSb) 

(7.6-11.3  /xm)  (PC  HgCdTe) 

2.0  X  10' W/cm^ 

1.0  X  10  W/cm^ 

Frame  rate 

2  forward  scans/sec 

Scanning  speed  during  data  acquisition 

36° /sec 

Detector  element  dwell  time 

520  /xsec 

Azimuth  position  resolution 

0.096  mrad 

Mirror  shaft  angle  encoder 

16  bits 

Total  incremental  count  (2.8°  in  azimuth) 

512  samples 

Dynamic  range 

12  bits/sample 

Sample  factor 

3.44  samples/dwell 

Internal  reference  cold 
hot 

Narcissus,  77K 
precision,  temp,  controlled 

|ChanneI  1  of  the  7.6-11.3  ;im  band  is  inactive. 

‘Internal  reference  source  occupies  -0.6°  of  the  azimuth  FOV,  leaving  -  2.2°. 
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Table  A-2  Summary  of  6MAP  data. 


No. 

Scene 

description 

Line-of-sight' 

Scan 

mode^ 

IRIG 

local 

time’ 

Frame 

(No.) 

Temp. 

(°F) 

RH 

(%) 

Visibility 

(mi) 

Wind 

Sun 

1 

Az 

El 

Speed 

Dir 

Az 

El 

Montauk  Point,  Aug  10,  1983 

1 

Overcast  sky 

227 

05 

8/8 

1-50 

75 

59 

10,000 

N 

110 

44 

feet/sec 

Bedford,  Sep  13,  1984 

1 

Hole  in  cloud  deck 

108 

12 

A 

3.8.52, 

1- 

66 

48 

>31 

calm 

_ 

112 

28 

21.3 

61 

2 

Hole  in  cloud  deck 

112 

15 

A 

3,8.55, 

62- 

66 

48 

>31 

calm 

112 

28 

36.8 

163 

3 

Hole  in  cloud  deck 

109 

19 

A 

3,8,58, 

164- 

66 

48 

>31 

calm 

_ 

113 

29 

39.8 

265 

4 

Hole  in  cloud  deck 

130 

11 

A 

3,9,07, 

266- 

66 

48 

>31 

calm 

115 

30 

51.3 

347 

5 

Hole  in  cloud  deck 

113 

09 

8/8 

3,9,11, 

348- 

66 

48 

>31 

calm 

_ 

116 

31 

47.8 

389 

6 

Hole  in  cloud  deck 

93 

08 

A 

3,9,18, 

390- 

67 

53 

>31 

calm 

_ 

117 

32 

57.8 

461 

7 

Hole  in  cloud  deck 

120 

14 

8/8 

3,9,38, 

521- 

67 

53 

>31 

calm 

121 

35 

38.8 

551 

'Az  =  degrees  clockwise  from  geographic  north;  El  =  degrees  above  horizontal. 


^Scan  mode:  A  =  alternate  mid-IR  and  long-lR  frames;  8/8  =  long-IR  channels  (5-12)  in  data  channels  (1-8) 
together  with  mid-lR  channels  (5-12)  in  data  channels  (9-16). 

^IRIG  time:  Day,  hour,  minute,  seconds. 
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Appendix  B 


AFGL  Data  Description 


The  information  provided  here  is  on  the  Air  Force 
Geophysics  Laboratory  (AFGL)  scenes.  These  scenes 
were  collected  from  an  IR  FLIR  Laboratory.  The 
Honeywell  IR  FLIR  is  intended  primarily  for  collect¬ 
ing  signature  data  but  has  also  been  applied  to  col¬ 
lecting  background  data.  The  principal  characteristics 
of  these  measurements  are  high  spatial  resolution, 
moderately  large  images,  and  a  wide  variety  of  avail¬ 
able  scenes.  A  description  of  the  instrument  is  given 
in  Table  B-1. 


Characteristics  of  AFGL  data  are  given  in  Table 
B-2.  All  scenes  contain  clouds  and  are  taken  with  fil¬ 
ter  #9  (7.84  to  13.16  ^m).  Table  B-3  gives  calibration 
data  for  all  scenes.  Note  that  scenes  designated  121-129 
have  a  smaller  FOV  from  typical  AFGL  FLIR  data. 
This  higher  resolution  is  achieved  by  using  an  objec¬ 
tive  lens  with  longer  focal  length  to  increase  spatial 
resolution.  The  system  spectral  response  is  given  in 
Fig.  B-1;  overall  response  is  approximately  80*l'’o  of 
that  achieved  without  a  filter  within  the  spectral 
bandpass. 


Table  B-1  Characteristics  of  the  AFGL-FLIR. 


Description 

Value 

Collecting  aperture 

25.2  cm^ 

Effective  focal  length 

7.67  cm 

System  //# 

1. 22 

Detectors 

Two  9-element,  serially  scanned 
PC  HgCdTe  arrays 

Detector  D*  (11.5  ^m,  lO  kHz) 

7.6  X  lO'"  cm  Hz’  VW 

Noise-equivalent  differential  temperature 

-0.2  K 

Spectral  bands 

lO-position  filter  wheel 

(APL  data  are  all  7.84-13. 16  ^m) 

Scanner 

10-facet  mirror  (horizontal), 
galvanometer  (vertical) 

Dynamic  range 

100:1 

Instantaneous  FOV 

0.43  mrad  (V)  x  0.36  mrad  (H) 

Total  FOV 

4.93°  (V)  X  7.62°  (H) 

Field  rate 

60  Hz 

Active  scan  lines 

192 

Video  output 

EIA  RS-170 

Table  B-2  Summary  of  AFGL  infrared  imagery  characteristics, 
(footnotes  for  this  table  follow  Table  B-3) 
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Table  B-3  Calibration  and  spectral  band  data  for  AFGL  infrared  imagery, 
(footnotes  for  this  table  and  Table  B-2  follow) 


Image 

Filter 

Spectral 

Band*’’ 

(/zm) 

Conversion  Factors*” 
to  Radiance 
(/iW/cmVsr) 

Image 

Multiplicative  Additive 

101 

19005A 

7.84-13.16 

8.4088- 

10“* 

1.8064 

10"’ 

101 

102 

19005A 

7.84-13.16 

5.3056- 

10* 

2.9287 

10"’ 

102 

103 

19005A 

7.84-13.16 

2.1122- 

10’ 

1.4909 

10"’ 

103 

104 

19005A 

7.84-13.16 

2.1122- 

10"’ 

1.4909 

10"’ 

104 

105 

19005A 

7.84-13.16 

8.4088- 

10* 

1.9456 

10"’ 

105 

106 

19005A 

7.84-13.16 

1.3327- 

10"’ 

1.4152 

10"’ 

106 

107 

19005A 

7.84-13.16 

1.3327- 

10"’ 

1.4126 

10"’ 

107 

108 

19005A 

7.84-13.16 

3.3476- 

10"’ 

-  2.4370 

10"’ 

108 

109 

19005A 

7.84-13.16 

1.0586- 

10"’ 

2.0466 

10"’ 

109 

110 

19005A 

7.84-13.16 

1.0586- 

10"’ 

2.5241 

10"’ 

110 

111 

19005A 

7.84-13.16 

8.4088- 

10“* 

1.0573 

10"’ 

111 

112 

19005A 

7.84-13.16 

8.4088- 

10"* 

1.3923 

10"’ 

112 

113 

19005A 

7.84-13.16 

1.6778- 

10"’ 

2.9420 

10"“ 

113 

114 

19005A 

7.84-13.16 

1.6778- 

10"’ 

4.1807 

10"“ 

114 

115 

B2180A 

7.96-11.39 

(4) 

(4) 

115 

116 

19005A 

7.84-13.16 

1.6778- 

10"’ 

-1.3316 

10"’ 

116 

117 

19005A 

7.84-13.16 

1.6778- 

10"’ 

5.6492 

10-“ 

117 

118 

19005A 

7.84-13.16 

1.6778- 

10"’ 

9.6727 

10-“ 

118 

119 

19005A 

7.84-13.16 

1.6778- 

10"’ 

9.6727 

10"“ 

119 

120 

19005A 

7.84-13.16 

8.4088- 

10* 

1.4649 

10-’ 

120 

121 

I9005A 

7.84-13.16 

6.6793- 

10"* 

3.2448 

10-’ 

121 

122 

I9005A 

7.84-13.16 

6.6793- 

10"* 

3.2414 

10-’ 

122 

123 

19005A 

7.84-13.16 

6.6793- 

10* 

3.2414 

10"’ 

123 

124 

19005A 

7.84-13.16 

8.4088- 

10* 

2.9607 

10"’ 

124 

125 

19005A 

7.84-13.16 

8.4088- 

10'* 

2.8066 

10"’ 

125 

126 

19005A 

7.84-13.16 

8.4088- 

10* 

3.2804 

10"’ 

126 

127 

19005A 

7.84-13.16 

8.4088- 

10"* 

3.2804 

10"’ 

127 

128 

19005A 

7.84-13.16 

8.4088- 

10* 

2.9684 

10"’ 

128 

129 

19005A 

7.84-13.16 

8.4088- 

10* 

3.1648 

10"’ 

129 

130 

19005A 

7.84-13.16 

1.0586- 

10"’ 

-  1.1037 

10"’ 

130 

131 

19005A 

7.84-13.16 

1.0586- 

10"’ 

-6.1401 

10“ 

131 

Footnotes  for  AFGL  background  data  description,  Tables  B2  and  B3 


'  FLIR  angular  resolution  and  digitized  sample  size  are  given  as 
value  in  the  horizontal  (scan)  direction  followed  by  value  in  ver¬ 
tical  direction.  The  images  are  formed  by  combining  the  (delayed) 
outputs  of  two  serially  scanning  detector  arrays. 

^The  digital  values  (0-2SS)  describing  the  brightness  of  each  pix¬ 
el  are  converted  to  actual  scene  radiance  (MW/cm^/sr/|<m)  by 
multiplying  the  digital  value  (DV)  of  the  data  by  the  multiplica¬ 
tive  factor  (MF)  and  adding  the  additive  factor  (AF): 


Brightness  =  MF  x  DV  +  AF. 

The  scene  has  a  maximum  of  64  gray  levels  at  multiples  of  4 
(i.e.,  0,  4,  8,...,  248,  252).  Therefore  the  quantization  level  is  4 
X  MF. 

’The  spectral  band  is  described  by  the  50*l7o  spectral  response 
points.  The  AFGL  FLIR  has  a  10-position  filter  wheel.  Refer¬ 
ence  1 1  gives  spectral  response  data  for  all  the  filters. 
'Calibration  data  supplied  for  image  115  are  not  correct. 
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Wavelength  (fim) 

Figure  B-1  Spectral  response  of  AFGL  FLIR  with  7.84-13.16^01  filter. 
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Appendix  C 


ERIM  Data  Description 


Information  on  the  Environmental  Research  Insti¬ 
tute  of  Michigan  (ERIM)  scenes  is  provided  in  this  ap¬ 
pendix. 


Data  provided  by  ERIM  consist  of  surface  scenes 
collected  by  the  M-7  multispectral  line  scanner  or  its 
similar  predecessor  the  M-S.  A  typical  collection  run 
is  made  at  low  altitude  from  a  C-47  aircraft.  The  scan¬ 
ner  looks  downward  (or  tilted  forward)  and  is  scanned 
(with  a  mirror)  in  the  cross-track  direction  over  a  90* 
sector.  Down-track  data  are  gathered  as  the  aircraft 
flies  over  the  scene.  The  principal  benefit  of  the  ERIM 
data  is  the  large  number  of  available  spectral  bands 
per  scene  and  the  large  size  of  some  scenes 
(>1,000,000  pixels). 


Table  C-1  gives  the  characteristics  of  the  M-7  scan¬ 
ner.  Table  C-2  summarizes  the  scene  data  provided 
to  APL.  Calibration  constants  to  convert  digital  counts 
to  radiance  or  apparent  temperature  (depending  on 
spectral  band)  are  provided  in  Table  C-3.  Spectral 
characteristics  of  the  major  detector  channels  are  given 
in  Fig.  C-1 .  There  are  two  infrared  detector  positions 
in  the  M-7  scanner.  One  is  devoted  to  the  long  wave¬ 
length  band  and  the  other  to  the  mid-IR  detectors. 
Several  mid-IR  detector  arrays  were  used  to  collect  the 
data.  Table  C-4  gives  the  relative  geometry  of  the  four 
detector  arrays  used  to  collect  these  data  and  data 
shifts  required  to  bring  pixels  of  different  spectral 
channels  into  approximate  alignment  in  the  center  of 
a  scan  line.  Note  that  geometric  effects  will  not  allow 
registration  across  an  entire  (single)  scan  line. 


Table  C-1  Characteristics  of  the  ERIM  M-7  scanner. 


Parameter 


Value 


Collecting  aperture 


5-inch 


Spatial  resolution  2  mrad 

Scan  rate  60/sec 


NEAT 


0.1  °C  nominal 


Spectral  bands  (^m) 


1.0- 1.4 
2.0-2.6 
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Footnotes  for  ERIM  background  data  description,  Tables  C-2  and  C-3 


I 

'  Angular  resolution  and  ground  sample  distance  given  as  value 
in  scan  (cross-track)  direction  followed  by  value  in  flight  (in-track) 

I  direction.  Resolutions  given  are  for  middle-infrared  bands;  spe¬ 

cifically  the  resolution  of  long-wavelength  bands  that  differ  from 
values  in  the  table  are; 


File 

Channel 

Spectral  band 

Resolution 

S2I 

4 

8.0-13.5 

3.5  mrad 

S23-S26 

4 

9.0-11.4 

2.9  mrad 

S29,F03, 

F05-F06 

3 

5.1-5.7 

2.9  mrad 

F04, 

F07-F09 

3 

9.0-11.4 

2.9  mrad 

^The  digital  values  (0-2SS)  describing  the  brightness  of  each  pix¬ 
el  is  converted  to  actual  scene  radiance  (/iW/cm^/sr/ftm)  or  ap¬ 
parent  temperature  (K)  by  multiplying  the  digital  value  (DV)  by 
the  multiplicative  factor  (MF)  and  adding  the  additive  factor  (AF): 

Brightness  =  MF  x  DV  -t  AF. 


The  units  of  the  brightness  result  are  given  in  the  units  column. 
MF  is  the  data  quantization  increment. 

Spatial  registration  is  given  as  the  relative  pixel  number  in  a  scan 
line  that  brings  the  images  of  different  channels  into  approxi¬ 
mate  spatial  alignment.  Alignment  is  a  function  of  the  scan  an¬ 
gle:  best  alignment  is  achieved  in  the  center  of  scan.  See  Table 
3,  Fig.  4,  and  Ref.  10  for  more  information. 

‘Channel  3  data  are  preferred  to  channel  2  data  for  the  9.0-11.4 
nm  spectral  band. 

Data  for  the  S. 1-5.7  jim  channel  are  noisy  and  may  not  be 
usable. 

Data  file  contains  discontinuous  images: 

F03  and  F04  contain  small  segments  (approximately  IS  scan 
lines  per  segment)  of  a  much  larger  image. 

F07  has  two  subimages  of  449  scan  lines  each. 

Tilted  scanner  images  have  vignetting  on  the  edges  of  the  image. 
Only  center  4(X)  pixels  are  completely  free  of  vignetting  effects. 
Original  channel  2  (3.0-4. 3  /.m)  of  image  F03  has  been  dropped 
(bad  data);  other  channels  were  renumbered.  Data  available  for 
1.0-1.4, 2.0-2.6,  and  8.0-13. S  ^m  channels  were  not  provided  to 
APL. 

'Long- wavelength  channels  of  M-S  scanner  cannot  be  easily 
brought  into  spatial  alignment  with  other  channels. 


Figure  C*1  Spectral  characteristics  of  ERIM  M-7  de¬ 
tector  channels. 
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Table  C-4  ERIM  middle-infrared  detector  arrays. 


Array*" 

Spectral  Bands*  ’ 

Images 

Detector  1 

Detector  2 

Detector  3 

InSb  3-3 

2.0-2.6  /im 

1.0-1. 4  /im 

4.5-5. 5  /im 

S21-S22 

InSb  3-7 

3. 9-4. 7  /tm 

2.0-2.6  /im*" 

3. 5-3. 9  /im 

S27-S30.  F04-F09 

InSb  3-8 

3.0-4. 2  /im 

2. 0-2. 6  /tm 

4.5-5.5  /im 

S01-S04,  S23-S26.  F03 

InSb  3-9 

3. 5-3.9  /im 

2.0-2.6  /tm 

4. 5-5. 5  /tm 

S05-S20 

'Detectors  are  circular  with  center-to-center  spacing  of  2.5  times  the  detector  diameter. 
^5. 1-5.7  /im  and  long-wavelength  detectors  are  at  a  different  focal  plane  of  the  M-7  scan¬ 
ner  and  are  in  spatial  registration  with  the  center  detector  of  the  above  arrays. 
^Detector  inoperative. 

Illustration  of  detector  shifting  to  achieve  spatial  coincidence 


Detector  footprints 
at  successive  samples 

Pixels  in  spatial  registration 
at  the  indicated  point 

r“ 

0  0  □ 

Detector  1 

(pixel  i  -  2) 

T 

0  0  0 

(pixel  i-  1) 

i 

m 

0  0  0  direction— 

Detector  2 

(pixel  i) 

e 

1 

0  0  0 

(pixel  i-f  1) 

0  0  0 

L_ 

Detector  3 

(pixel  i  -t-  2) 
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Figur*  D-S  Power  spectral  densities  of  ERIM  scenes 
of  Camp  A.P.  Hill,  8.0-14.0  ^m  spectral  band,  (a)  Morn¬ 
ing  (S01CH3).  (b)  Afternoon  (S02CH3).  (c)  Evening 
(S03CH3). 


Spatial  frequency  (cycles/mrad) 

Figure  D-6  Power  spectral  densities  of  ERIM  scenes 
of  Port  Hueneme  (port  and  water),  (a)  4.5-5.5  spec¬ 
tral  band  (S23CH3).  (b)  9.0-11.4  pm  spectral  band 
(S23CH4). 
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I  10-9 


I  10-9 


^o-i2i _ i..miua _ I  I  _ I  I  mini 

io-«  10-3  10-2  10-’ 


Spatial  frequency  (cycles/mradl 


Figur*  E-4  Power  spectral  density  of  ERIM  scene  of 
Port  Hueneme,  spectral  band  4.5-5.5  after  high- 
pass  filter  (scene  S23CH3). 


io-’2l_ 

10-<  10-3  10^2  10-1  10° 

Spatial  frequency  (cycles/mrad) 

Figure  E-3  Power  spectral  densities  of  ERIM  after¬ 
noon  scene  of  Camp  A.P.  Hill  after  high-pass  filter,  (a) 
2.0-2.6  pm  spectral  band  (S02CH1).  (b)  4.5-5.5  spec¬ 
tral  band  (S02CH2).  (c)  8.0-14.0  /im  spectral  band 
(S02CH3). 
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